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Abstract 

This paper concerns the existence of global weak solutions to the barotropic 
compressible Navier-Stokes equations with degenerate viscosity coefficients. We 
construct suitable approximate system which has smooth solutions satisfying the 
energy inequality, the BD entropy one, and the Mellet-Vasseur type estimate. 

Then, after adapting the compactness results due to Bresch-Desjardins (2002, 2003) 
and Mellet-Vasseur (2007), we obtain the global existence of weak solutions to the 
barotropic compressible Navier-Stokes equations with degenerate viscosity coeffi¬ 
cients in two or three dimensional periodic domains or whole space for large initial 
data. This, in particular, solved an open problem proposed by Lions (1998). 
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1 Introduction and main results 

The barotropic compressible Navier-Stokes equations, which 
scribing the evolution of a viscous compressible fluid, read as 

p t + div(ptt) = 0, 

(pu)t + di v(pu <g) u) — div§ + VP(p) = 

where x € U C M N (N = 2,3 ),t > 0, p is the density, u = (u\, ■ ■ ■ ,un) is the velocity, 
§ is the viscous stress tensor, and P(p) = ap' y (a > 0,7 > 1) is the pressure. Without 
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are the basic models de- 
follows 
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( 1 . 1 ) 
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loss of generality, it is assumed that a = 1. Two major cases will be considered: either 


§ = §1 = hVu + gdivitl, 

( 1 . 2 ) 

S = §2 — hVu + gdivul, 

(1.3) 


where Vu = ^(Vu + (Vrt) tr ), I is the identical matrix, and h,g satisfy the physical 
restrictions 

h > 0, h + Ng > 0. (1.4) 

There are many studies on the global existence and behavior of solutions to (11.111 
(11.21) when both h and g are constants. The one-dimensional problem has been studied 
extensively, see |15ll27ll28] and the references therein. For the multi-dimensional case, 
the global classical solutions with the density strictly away from vacuum were first 
obtained by Matsumura-Nishida [34] for initial data close to a non-vacuum equilibrium. 
Recently, Huang-Li-Xin m obtained the global classical solutions with the density 
containing vacuum provided the initial energy is suitably small. For the weak solutions, 
Hoff jl 6 Hl 8 ] studied the problem for discontinuous initial data. When the initial total 
energy is finite (which implies that the initial density may vanish), Lions j31] obtained 
the global existence of weak solutions provided the exponent 7 is suitably large, which 
was further relaxed by Feireisl-Novotny-Petzeltova m to 7 > 3/2 for three-dimensional 
case. 

On the other hand, there are important and interesting phenomena where h and g 
depend on the density which are degenerate at vacuum. Indeed, as pointed out by Liu- 
Xin-Yang in [32], in the derivation of the compressible Navier-Stokes equations from 
the Boltzmann equation by the Chapman-Enskog expansions, the viscosity depends 
on the temperature, which is translated into the dependence of the viscosity on the 
density for barotropic flows. Moreover, Lions m also proposed various models for 
shallow water, in particular, he points out that the global existence of weak solutions 
to (11.11) (11.21) with h = p, g = 0 remains open. Recently, a friction shallow-water 
system, with flat bottom topography, which is derived in [8[[9[[T2[[33] . can be written 
in a two-dimensional space domain H as EH) E3D with h(p) = g(p) = p. Indeed, 
such models appear naturally and often in geophysical flows [EEIEHS]- Therefore, it is 
of great importance to study the compressible Navier-Stokes equations ( 11 . 11 ) ( 11 . 21 ) and 
(11.11) (11.31) with density-dependent viscosity. 

In the one-dimensional case with h = g = Ap a for some positive constants A and 
a, the well-posedness of either the initial value problem or the initial boundary value 
ones with fixed or free boundaries has been studied by many authors (see [ 22il30 L32l 
[36l[39llTTI443j and the references therein). In higher dimensions, assuming that h is a 
constant and g(p) = ap° with a > 0 and (3 > 3, Vaigant-Kazhikhov [37] first proved 
that for the two-dimensional case EH (ESI) with slip boundary conditions has a unique 
global strong and classical solution. Recently, for the Cauchy problem and the periodic 
boundary conditions, Huang-Li [191120] and Jiu-Wang-Xin [241426] relaxed the condition 
f3 > 3 to f3 > 4/3. For the case h = h(p) and g = g(p ), in addition to (11.41) . under the 
condition that 

g(p) = ti(p)p-h(p), (1.5) 

Bresch-Desjardins [3HZ] have made important progress. Indeed, for the periodic bound¬ 
ary conditions and the Cauchy problem, they succeeded in obtaining a new entropy 
inequality (called BD entropy) which can not only be applied to the vacuum case but 
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also be used to get the global existence of weak solutions to (EU) (E 2 ) and CCD OK 
with some additional drag terms m\- Later, by obtaining a new apriori estimate on 
smooth approximation solutions, Mellet-Vasseur [35] study the stability of barotropic 
compressible Navier-Stokes equations ( 11 . 11 ) ( 11 . 21 ) and (11.11) (11.31) without any additional 
drag term. However, the construction of the smooth approximation solutions remains 
to be carried out, which does not seem routine in the case of appearance of vacuum. In 
fact, only part results for special cases are available. In particular, for one-dimensional 
case, Li-Li-Xin m obtained the global existence of weak solutions to (11.11) (II .31) with 
h(p) = g(p) = p a (a > 1 / 2 ) and proved that for any global entropy weak solution, 
any vacuum state must vanish within the finite time. Later, when the initial data is 
spherically symmetric, Guo-Jiu-Xin m obtained the global existence of weak solutions 
to (II. ip (II.2p whose Lagrange structure and dynamics are studied by Guo-Li-Xin [T31j . 
Thus, the main aim of this paper is to obtain the global existence of weak solutions to 
(1LTD (fL2D and (jl.ll) (jl.3l) for 7 > 1 and for general initial data by constructing some 
suitable smooth approximation solutions. 

For the sake of simplicity, it is assumed that for constant a > 0, 

Kp) = P a , 9{p) = (a - 1 )p a . ( 1 - 6 ) 

We then consider the Cauchy problem, f l = R n (N = 2,3), and the case of bounded 
domains with periodic boundary conditions, H = T n (N = 2,3). The initial conditions 
are imposed as 

p(x, t = 0) = p 0 , pu(x, t = 0) = m 0 . (1.7) 

We always assume that the initial data po,mo satisfy that for some constant r/o > 0, 

po > 0 a.e. in H, po ^ 0, po € L 1 (fl) n L 7 (H), Vpg -1 ^ 2 G L 2 (H), 

< mo € L 27 A 7+1 )(H), 777-0 = 0 a.e. on Ho, (1-8) 

^Po^KI 2 ^ 0 el 1 ^), 

where we agree that p 0 1 Vo \mo\ 2+r10 = 0 a.e. on Ho, the vacuum set of po, defined by 

Ho = {x € H \po(x) = 0}. (1.9) 

Before stating the main results, we give the definition of a weak solution to m 
(11.31) (11.61) (11.71) . Similarly, one can define a weak solution to (11.11) (11.21) (11.61) (11.71) . 

Definition 1.1 For N = 2,3, let H = T N or H = (p, u) is said to be a weak 

solution to (11.11) () 1.3 D (11.61) (11.71) if 

'0 < p G L°°(0, T ; L X (H) n L 7 (H)), 

V p( 7 +a-i)/2 € l 2 (0,T; (^(H))^), 

< Vp“ -1 / 2 , ^fpu G L°°(0, T; (L 2 (H)) Ar ), 
h(p)Vu, /7(p)(Vn) tr €L 2 (0,T-,(W^'\n)) NxN ), 

,9(p) divu G L 2 (0,T; W^’^H)), 

with (p, yfpu) satisfying 

f« + div(^») = 0, mV , (1.10) 

\p(x,t = 0) = p 0 (x), 
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and if the following equality holds for all smooth test function <f(x, t ) with compact 
support such that f>(x, T) = 0 : 


/ mo • 4>(x, 0)dx + / / + y/pu <g) yfpu : V0 4- p^divcf) dxdt 

Jn Jo Jn 

~ \{Kp)^ u ^ ( 1 ) ) - u ) tr ■> V<£) - (fl'(p)divu, div</>) = 0, 


( 1 . 11 ) 


where 


(h(p)Vu, V0) = - f [ p a ~ l/2 yfpu • 

Jo Jn 

I f y/pujdip a - l / 2 di(t)jdxdt , 
Jo JO 


2a - 1 

rT 


(h(p)(f\7 u) tr <p) = — f f p a x ^ 2 yfpu ■ Vdivcfdxdt 

Jo Jn 

sjpuidj p a ~ l ^ 2 di <fj dxdt , 


2a 

2a - 1 


'o Jn 


(g(p)divu, div</>) = — (a — 1) f f p a x ^ 2 yfpu ■ Vdivfdxdt 

Jo Jn 

— — -11 f f . Vp Q_1/ 2 divcfdxdt. 

2a — 1 J 0 J n 

Then the first main result of this paper is as follows: 

Theorem 1.1 Let fl = R 2 or T 2 . Suppose that a and 7 satisfy 

a > 1/2, 7 > 1, 7 > 2a — 1. 


( 1 . 12 ) 


Moreover, assume that the initial data (po,mo) satisfy (11.81) . Then there exists a global 
weak solution ( p,u ) to the problem (11.11) (11.31) ([1.611 (jl.7|) . 


The method of Theorem o can be applied directly to the system (HI]) HI]), that 


is 


Theorem 1.2 Let Ll = R 2 or T 2 . Under the conditions of Theorem M.il there exists a 
global weak solution (p, u ) to the problem ( 11 . 11 ) ( 11 . 21 ) () 1 . 6 1 ) ( 11 . 71 ) . 

Theorems 11.11 and OI are concerning with the two-dimensional case. As for the 
three-dimensional case, we have 

Theorem 1.3 Let Q = M 3 or T 3 . Suppose that a E [3/4,2) and 7 € (1,3) satisfy 


7 E 


(1,6a — 3), for a E [3/4,1], 

[2a — 1, 3a — 1], for a E (1, 2). 


(1.13) 


Assume that the initial data (po,mo) satisfy (11.81) . Moreover, if a E (1,2), in addition 
to (11.81) . we assume that 

Po 3 |m 0 | 4 E L 1 ^), (1.14) 

where we agree that /)Q 3 |mo| 4 = 0 a.e. on f2o as in (11.91) . Then there exists a global 
weak solution ( p,u ) to the problem (II. ip (11.21) ([1.61) (jl.7|) . 
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Theorem 1.4 Let 12 = M 3 or T 3 . Suppose that a = 1 and 7 E (1,3). Assume that the 
initial data ( po,mo ) satisfy (11.81) . Then there exists a global weak solution ( p,u) to the 
problem (11.11) (ll.3p (11.61) (|1.71) . 

A few remarks are in order: 

Remark 1.1 If a = 1 and 7 = 2 , Theorems 0[Z2 give a positive answer to the 
open problem proposed by Lions JH Section 8 . 4 ]: “In the first case (ie (8.70)-(8.71)), 
the Cauchy problem is completely open for the models involving (8.73)”, where (8.70)- 
(8.71) (8.73) is corresponding to (11.11) (11.21) (11.61) with a = 1 and 7 = 2 . 

Remark 1.2 For three-dimensional case, it should be noted that Theorem P here is 
valid for all 7 € (1,3) provided h = p and g = 0. Therefore, for h = p and g = 0, our 
Theorems ro and \l-4\ establish the existence of global weak solutions to mu m and 
(ED El with 7 E (1,3) for general initial data, which is in sharp contrast to the case 
that h and g are both constants, where the condition 7 > 3/2 is essential in the analysis 
of Lions m and Feireisl-Novotny-Petzeltova m ■ In fact, for h and g being both 
constants and 7 E (1,3/2], it remains completely open to obtain the global existence of 
weak solutions to (O) dLZD for general initial data except for the spherically symmetric 
case JM). 

Remark 1.3 After some routine modifications, for the system (11.11) (11.21) . our method 
can be applied directly to the case that h and g satisfy (11.51) and the conditions listed in 
m together with some additional constraints. However, for the shallow-water system 
with h = g = p ( ;m. 12\\3ffi). the global existence of weak solutions to ED EH) or 
dnp eh for general initial data remains open since (USD fails for this case. 

Remark 1.4 For the system (ED El and three-dimensional case, our construction 
depends on the condition that h = p and g = 0, and cannot be applied directly to the 
general case that h and g satisfy the conditions listed in Theorem 11.31 This will be left 
for future. 

Remark 1.5 Around the same time when this paper is finished, there are announce¬ 
ments of some existence results on the problem EH) El El ED with a = l and 
12 = T N (N = 2,3) by Vasseur-Yu !38 l with a different approach. However, we have 
difficulties to understand some of their key a priori assumptions near vacuum in their 
arguments. 

We now make some comments on the analysis of this paper. Since the compactness 
arguments are similar to those of Bresch-Desjardins m and Mellet-Vasseur [35], the 
main point of this paper is to construct smooth approximate solutions, whose densities 
are bounded from above and strictly bounded away from vacuum provided the smooth 
initial ones are, satisfying the energy estimate, the BD entropy inequality, and the 
Mellet-Vasseur type estimate. To this end, we first deal with the periodic case and 
consider the following approximate system 

p t + di v{pu) = £p 1 / 2 div(p~ 1 / 2 h' e (p)'V p), 

< put + pu • Vu - dw(h £ (p)Vu) - \7(g £ {p)divu) + VP (1-15) 

„= y/Ed\v(h e (p)Y7u) + \/iV (g £ (p)divu) - e~ £ 3 (p £ ~ 2 + p~ e ~ 2 )u , 
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where 

Mp) = P° + e 1/3 (p 7/8 + P 7 ), 9e(p) = pK(p) ~ Mp), (1-16) 

with 

0 < e < e 0 — min{(2a — l)(16(a+ 7)) _10 ,?7o}, 7 = 7 + 1/6. (1-17) 

Here, we propose to approximate (II. lib by (11.1511 1 which is a parabolic equation for 
any fixed e > 0 and hence has smooth effects on the density provided the smooth 
initial density is strictly away from vacuum. The specific choices of the higher order 
regularization in (11.151) have several key advantages. First, it can be shown that the 
smooth solutions to the new system (11.1511 — (11.171) satisfy the energy and the Mellet- 
Vasseur type estimates. Moreover, after some careful calculations, we find that the 
most difficult term induced by sp 1 ^ 2 dw(p~ 1 ^ 2 h' e (p)'V p) has the right sign (see (12.171) 1 
which implies that the solutions to our approximate system also satisfy the BD entropy 
inequality. In fact, this is one of the key observations of this paper. Next, in order 
to obtain the lower and upper bounds of the density, in addition to the estimate on 
L°°(0, T; L Ar+<s )-norm of p x ^ N+ ^u which can be obtained for the system (11.11) (11.31) 
in two-dimensional case (see (I2.29j) l and for (jl.ll) (jl.2j) in both two-dimensional (see 
(12.291) 1 and three-dimensional cases (see (13.71) 1. one still needs some additional estimates 
on the L°°(0, T; L p )-norm (for suitably large p) of p and p -1 which can be achieved 
by adding a damping term — e~ £ ( p £ + p~ £ )u on the righthand side of (1 1.1 1) o (see 
(11.151) 9). However, for e —>• 0 + , this term will bring new difficulties which can be 
overcome by adding e 1//3 (p 7 / 8 + p 7+1 / 6 ) to h(p) (see (11.161) 1. This idea is motivated by 
our previous study on the one-dimensional problem m- With all these estimates at 
hand, we can use a De Giorgi-type procedure to bound the density from above and 
below, in particular, the density is strictly away from vacuum provided the initial one 
is (see (12.251) 1. In fact, this is another key issue of this paper. Once we obtained (|2.25|) . 
we can use the L p -theory for parabolic system to get the estimates on the L p (0, T; L p )- 
norm of (p,u), (pt,ut), and (V 2 p, V 2 u) (see (12.401) 1. This in turn implies that the 
approximate system (11.151) - (11.171) has a global strong solution with smooth initial data. 
Next, after adapting the compactness results due to Bresch-Desjardins m and Mellet- 
Vasseur [35], we can obtain the global existence of the weak solutions to either (11.11) 
(TOD (fL6l) for two-dimensional periodic case or ([LTD (fL2l) (fL6l) for two-dimensional 
and three-dimensional periodic cases. Finally, to prove Theorems 11.11 and 11.31 where 
n = R n (N = 2,3), some extra care should be taken due to the unboundedness of 
the domain. In this case, we consider the system (|1.15jl — (|1.17D in Q E = (—e _CT °, e~ a °) 2 
and (13.21) (11.161) (I1.17|) in Q e = (— e~ a °, e~ r7 °) N (N = 2,3) and impose the Neumann 
boundary condition on p and Navier-slip conditions on u (see (] 5.6 D and (|5.15l) 'l. Then 
we can adapt the preceding proofs in the case H = T n (N = 2,3) to H = R n (N = 2,3). 

This paper is organized as follows. Since the proof of Theorem 11.21 is similar as that 
of Theorem 1 1.1 1 after some routine modifications, we will only prove Theorems 11.11 fL3l 
and 11.41 In the next section, we work on the problem (11.11) (11.31) (11.61) (11.71) in the 
two-dimensional periodic case, H = T 2 , then in the Section 3, we adapt the previous 
procedure to the problem (11.11) (11.21) (11.61) (11.71) in the three-dimensional periodic case, 
n = T 3 . Next, in the section 4, we will construct a new approximate system which can 
be applied to obtain the global weak solutions to the problem (|l.ll) (jl.31) (11.61) (jl.7D in 
the three-dimensional periodic case, and in the end (Section 5) we shall explain how to 
modify the preceding proofs in the cases where H = R N (N = 2, 3). 
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2 Proof of Theorem 11.11 : 0, = T 2 


In this section, we study the 2-dimensional periodic case which is the simplest one, yet 
the most important case since most of ideas to be developed here can be modified to 
deal with other cases. 


2.1 A priori estimates 

For e as in (11.1711 , let smooth functions po £ > 0 and uq £ satisfy 

l|POe||l,inLT'(fi) + \\V Poe llz, 2 (n) + ^ ' ’ 11 ^ Po^ \\l 2 (O) + £l/ " II ^ Poe ^ ’11^(0) 

+ e 13 / 3 e £ ||Poe Vw+e 137 ^ E II Pq £ ^llz^n) ^ 2 - 1 ) 

+ e 4 e- £ - 3 ||p £ 0 ; 2+ “- 1 || il(n) + £ 4 e- £ - 3 ||p 0 -/- 2+Q - 1 |hi W < C, 

and 

f P 0 e\y- 0 e\ 2+VO dx < C, (2.2) 

Ju 

for some constant C independent of e. We extend po £ and u 0 e ^-periodically to M 2 and 
consider the system fll.15p - fll.17D with initial data: 

(p,u)(x,0) = (p 0e ,u 0e ). (2.3) 

Let T > 0 be a fixed time and (p, it) be a smooth solution to (ll.15jl - fjl.17jl (j2.3H on 

n x (o, t]. 

Then, we will establish some necessary a priori bounds for ( p,u ). The first one is 
the energy-type inequality. 


Lemma 2.1 There exists some generic constant C independent of e and T such that 

rT 
/ o 


sup [ {p\u\ 2 + p + p 1 )dx + e f fp l h' e {p)\'\7 p\ 2 {\+ \u\ 2 )dxdt 
<t<T J Jo J 

h £ (p)\Vu\ 2 dxdt + e~ £ 3 J J (p £ ” + p~ £ 2 )\u\ 2 dxdt < C, 


o <t<T 
rT 


where and throughout this section, for any f, 


! fdx -.L 


fdx = I fdx. 

n 

Proof. First, integrating (jl,15D 1 over 12 x (0, T) together with (12.31) gives 

rT r 

i r~ lu ' 

/ 0 


sup / pdx + £ / 
0<t<T J Jo 


P~%(p)\Vp\ 2 dxdt<C. 


(2.4) 


(2.5) 
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Next, multiplying (I1.15|) 2 by u, integrating by parts, and using (I1.15D 1 yield 
J p\u\ 2 dx)t + J(h £ (p)(\Vu\ 2 + y /£\Vu\ 2 ) + (l + y/ r e)g £ (p)(dwu) 2 )dx 

J p 1 ^ 2 dw(p~ 1 ^ 2 h' £ (p)'Vp)\u\ 2 dx 


= —- / p 1 h! £ (p)\S7p\ 2 \u\ 2 dx — e / h' £ (p)Vp ■ V u ■ udx 


( 2 . 6 ) 


< - 


< - 


^ J p 1 h' £ (p)\Vp\ 2 \u\ 2 dx + 2e J ph' £ (p)\V u\ 2 dx 

| J P -1 ^(p)|Vp| 2 M 2 dx + ^ J h £ (p)\Vu\ 2 dx. 


Then, to estimate the last term on the left hand side of (I2.6p . after integration by 
parts and using (|1.15D 1 . one obtains that for q ^ 1, 


u ■ V p q dx = — 


q 


p g Mi v(pu)dx 

p q ~ l (-pt + £p l/2 d\v(p~ 1/2 h' £ (p)V p))dx (2.7) 


<7-1.7 

ihj 

1 ' r P9dx)t+ 9 2(q-l) J P^KiP^Pl 2 ^- 


q -1 

Finally, for u € M Ar (iV = 2, 3), we have 

(divu) 2 < N\Vv\ 2 < lV|Vu| 2 , 

which together with () 1.16 f) implies that for N = 2,3, 

f d(h £ (p)\Vv\ 2 + g e {p){ divu) 2 ) > min{ Na - (N - 1), 1 }h £ (p)\Vv\ 2 , 
[ 4(h e (p)|Vu| 2 + g £ {p)( divu) 2 ) > minjlVa — (N — 1), l}h £ (p)\Vv\ 2 . 


( 2 . 8 ) 


Since e < eq, the combination of (12.51) - (12.811 . (12.11) . with (12.21) yields (12.41) . which com¬ 
pletes the proof of Lemma 12.11 

Now we are in a position to derive the following entropy estimate which in particular 
yields the uniform BD one due to Bresch-Desjardins m\- 

Lemma 2.2 There exists some generic constant C independent of e and T such that 

sup [ (p-\h' e (p)) 2 Np\ 2 + e 13/3 e~ £ ~ 3 (p 6 ' 2 ^- 1 + p-^ 2 - 1 / 8 )) dx 
KKTJ v ' 

h £ (p)\Vu\ 2 dxdt + J J p y ~ 3 h e (p)\Vp\ 2 dxdt < C. 


0<t<T . 

+ f 


(2.9) 


Proof. First, set 


G^ £/ 9 1 / 2 div(p- 1 / 2 h'(p)Vp) 


( 2 . 10 ) 


and 


p's (p) = p-%(p) > 0. 


(2.11) 














Multiplying (11.151) i by (p £ (p) leads to 

(<Pe(p))t + u ■ Vp £ (p) + p(p £ (p)divu = p' e (p)G , 

which gives 

(S7p e (p))t + u ■ VVp £ (p) + Vu • V(p e (p) + V {pp' £ (p)dWu) = V(ip £ (p)G). 
Thus, multiplying (12.121) by pVp £ (p) and integration by parts show that 


( 2 . 12 ) 


1 


(/ p\X7p £ (p)\ 2 dx)t + / Vh £ (p) -Vu -Vp £ (p)dx 


+ / Vh £ {p) ■ V(p<p £ (p)divu)dx 


(2.13) 


+ 


J v'e(p)G h £ (p) - \p' £ {p)\^p\ 2 ^j dx = 0. 

Next, multiplying (11.151) 9 by S/<p e (p) leads to 

J UfVh £ (p)dx + J u ■ Vu ■ Vh £ (p)dx — (1 + y/s) j h £ (p)Vdxvu ■ V tp £ (p)dx 
- (1 + y/e) J Vh £ (p) ■ Vu • Vp £ {p)dx + (1 + y/e) J g £ {p)ddvuAip E {p)dx 

+ j p \pWAp)\Vp\ 2 dx + e~ £ 3 J ( p £ 2 + p~ £ 2 )u ■ Vp £ (p)dx = 0, 
where the following simple fact has been used: 

— J h £ (p)Au ■ Vp E (p)dx = — J h £ (p)V divu ■ Vp £ (p)dx. 

Since (11. 15 D i implies 

{h e (p))t + di v(h £ (p)u) + {ph' e (p) - h £ {p))divu = h £ (p)G, 
the first term on the left hand side of ()2.141) is handled as 

J u t ■ Vh £ (p)dx = (J u- Vh £ (p)dx) t - J u ■ Vh £ (p) t dx 

= {J u ■ Vh £ {p)dx) t - / u ■ Vu ■ Vh £ (p)dx 
— 2 j h £ (p)Vu : Vudx + J h £ (p)\Vu\ 2 dx 
— J (ph £ (p) — h £ (p))(divu) 2 dx + J dwuh £ (p)Gdx, 
where in the second equality one has used 
J u ■ V d\v{h £ {p)u)dx 
= - J diU ■ V{h £ (p)ui)dx 
= — J u ■ Vu ■ Vh £ (p)dx — J h £ (p)diU ■ V'Uidx 
= — J u ■ Vu ■ Vh £ (p)dx — 2 J h £ {p)Vu :Vudx + J h £ {p)\Vu\ 2 dx. 


(2.14) 


(2.15) 
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J 9e (p)divuAip £ {p)dx 

= J Vh £ (p) ■ V (pip' £ (p))divudx + J pip £ (p)Vh £ (p) ■ Vdivudx 

— J h £ (p)Vdivu ■ V(p £ (p)dx — J g £ {p)Vddvu ■ Vip £ (p)dx 

- J di vuVg e (p) ■ Vp £ (p)dx 

= J (V/i e (p) ■ V(py' £ {p)) - Vg e {p) ■ Vip £ (p)) di vudx 

+ J (w4(p)VMp) - h £ (p)V(p e (p) - g e (p)Vip £ (p)) • Vdi vudx = 0 

due to (11.161) and (12.111) . 

Since (12.111) and (12.101) imply 


(2.16) 


Now, multiplying (12.131) by 1 + y/e and adding the resulting equality to (12.141) . one 
can obtain after using (12.151) that 

P\VVe{p)\ 2 dx) t + (J pu-V<p £ (p)dx) t + J h £ (p)\V u\ 2 dx 

+ j p '(p)<p' £ (p)\^P\ 2dx + e ~ £ 3 J ( P £ 2 + P~ £ 2 )u-\7tp £ (p)dx 

+ (1 + Ve) J v' e {p)G h £ {p) - )pt e (p)\Vp\ 2 + x + pdivu'j dx 

= 2 j h £ (p)Vu : Vudx + J ( ph' £ (p ) — h £ {p))(ddvu) 2 dx 
<\j h £ (p)\V u\ 2 dx + C J h £ (p)\Vu\ 2 dx, 
where in the first equality one has used the following simple calculations: 

J Vh £ (p) ■ V(p(p' £ (p)divu)dx — J h £ (p)Vdivu ■ Vip £ (p)dx 
+ 


Ah e (p)--< f / e (p)\Vp\ 2 = £- 1 G, 


the last term on the left hand side of (|2.16l) satisfies 

(1 + y/e) J <-p's{p) g - ^' £ (p)\Vp\ 2 + x + ^- pdivuj dx 


2 dx — % 


p 2 tp' £ (p) (dmt)" dx 


> 


h J^ p)a 

— J (p' £ (p)G 2 dx — Ce J h e (p)\Vu\ 2 dx. 


(2.17) 
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Finally, it follows from (12.111) and (11.161) that 

_ 2 —2 _ 

(. P £ + P~ £ )u-Vp £ (p)dx 

-M? 

which, together with (12.71) . (12.161) . (12.171) . (12.41) . and (12.11) . yields (12.91) . The proof of 
Lemma 12.21 is finished. 

With Lemmas no and 12.21 at hand, we can prove the following Mellet-Vasseur type 
estimate ( [35]). 

Lemma 2.3 Assume that 7 > 1 satisfies 7 > (1 + a)/2 in addition. Then there exists 
some generic constant C depending on T but independent of e such that 

sup / p{e + \u\ 2 ) ln(e + \u\ 2 )dx < C. (2.18) 

0 <t<T J 


,e~ 2 +a -1 7 e 7/3 e- 2 -l /8 

H-5-1" 


7e 7//3 p e 2+ ^ 1 


2 + a — 1 

a / o _e_2+a_1 
—e -2 + a — 1 


— e* 1 + (7 — l)e 2 

7e 7 / 3 p- £ " 2 - 1 /8 7£ 7 / 3 p - £ - 2 +t-i' 


dx 


+ £ 2 


1 — (7 — l)e 2 


dx, 


Proof. First, multiplying (H.15l) o by (1 + ln(e + |«| 2 ))u and integrating lead to 
\~lt [ + l u l 2 ) l R ( e + \ u \ 2 )dx — i f (e + |«| 2 ) ln(e + \u\ 2 )Gdx 


+ 1 + ln(e + \u\ 2 ))(h £ (p)(\Vu\ 2 + v^lVul 2 ) + (1 + y/e)g e (p)(divu) 2 )dx 


<C h £ (p)\\7u\ 2 dx — / (1 + ln(e + |tt| 2 ))rt • Vp 7 dx 


< C 


J h £ (p)\\7u\ 2 dx + C J ln 2 (e + |ri| 2 )p 27 "dx, 


2 UiUk 


diUkfi^dx 


where in the last inequality one has used the following estimate 
(1 + ln(e + \u\ 2 ))u ■ S/p^dx 

< f (1 + ln(e + |'u| 2 ))|divii|p 7 dx + f „ 

J J e+\u\ z 

< C J \r?{e + \u\ 2 )p 2l ~ a dx + C J h £ (p)\Vu\ 2 dx. 

Then, integration by parts gives 

— i J(e + |u| 2 ) ln(e + \u\ 2 )Gdx 

= J fp-'h'M |Vp| 2 (e + !“| 2 ) ta(e + + 6 / K(P)V P ■ Vu • u dl 


> 


+ £ 
£ 


j h' £ (p ) ln(e + |«| 2 )Vp • Vtt • rtdx 
J p~ l h' e {p)\ Vp| 2 (e + |u| 2 ) ln(e + \u\ 2 )dx — £ J ph' £ (p)\Vu\ 2 dx 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


— £ / p l h £ (p)\Vp\ 2 \u\ 2 dx — 2e(a + 2) / h e (p) ln(e + |ri| 2 )|V , u| 2 dx. 
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It follows from this, (12.191) . (12.41) . (12.91) . (12.81) . and (12.21) that 


sup f p(e + \u\ 2 ) ln(e + \u\ 2 )dx < C + f f ln 2 (e + \u\ 2 )p 21 a dxdt. ( 2 . 22 ) 
5<t<T J Jo J 


o <t<T J JO 

Finally, since 7 > (a + l)/2, it holds that 


ln 2 (e + \u\ 2 )p 2 ' r a dx < C J(p + p 2l ){ 1 + \u\)dx 

<C + C J p\u\ 2 dx + C J (p + p 4 ' y ~ 1 )dx 


(2.23) 


<C, 


where in the last inequality, one has used (|2.4D . (12.91) . and the following Sobolev in¬ 
equality that for any p > 1 , there exists some constant C depending only on a and p 
such, that 

IWIlpPD < C\\p\\ LHSl) +C\\Vp a - 1 / 2 \\%^ ) ~ 1) . (2.24) 

Putting (12.231) into (|2.22l) yields (12.181) . The proof of Lemma 12431 is completed. 

Next, we will use a De Giorgi-type procedure to obtain the following estimates on the 
lower and upper bounds of the density which are the key to obtain the global existence 
of strong solutions to the problem (II . 15|) — (II. 1 7D (12.31) . 


Lemma 2.4 There exists some positive constant C depending on £ and T such that 
for all (. x,t ) € 12 X (0,T) 

C- 1 < p(x,t) <C. (2.25) 


Proof. First, multiplying (I1.15IL by \u\ £ u and integrating in space give 
1 d 


J p\u\ 2+£ dx — 2 ^ J \u\ 2+£ Gdx 


2 “b £ dt 

+ J \u\ £ (h £ (p)(\Vu\ 2 + y/s\\7u\ 2 ) + (1 + ^/s)g e (p)(divu) 2 ) dx 
£■(1 + 2-^/e) 


+ 


J h e (p) | 


u\ £ \\7\u\fdx + e 


(p £ + p £ )\u\ 2+£ dx 


= —- J h £ (p)u ■ Vu • X7\u\ £ dx — (1 + y/e) J g £ (p)divuu ■ X7\u\ £ dx 
— J \u\ £ u • Vp 7 dx 

< 4(a + 7 ) 5 - J h £ (p)\u\ £ \Vu\ 2 dx + C J h £ (p)\Vu\ 2 dx 
+ C f [p £ 2 + p~ £ \u\ 2 dx + C, 


where in the last inequality one has used the following simple fact that 

V" 2 + P~ 


sup 
0 <t<T 


' £ 2 1 -- £ 2 ' dx < C, 


(2.26) 


(2.27) 
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due to (12.91) . Integration by parts yields that 

- ih / M2+ ‘ Gdx 

p~ l h' e {p)\W p\ 2 \u\ 2+s dx + e J h' e (p)\u\ 1+£ 'Vp ■ V\u\dx 
■ J p~ 1 h' £ (p)\'V p\ 2 \u\ 2+£ dx — 2(2 + s)e J ph' e (p)\u\ e \\7u\ 2 dx. 


2(2 + s) 


~ 8(2 + e) 

It follows from (12.26j) . (j2.28j) . (12.81) . (12.41) . (12.9p . and (11.17[) that 


(2.28) 


sup / p\u\ 2+£ dx + yfe / / h £ (p)\u\ £ \Vu\ 2 dxdt < C. 

o <t<T J Jo J 

Next, since v = p 1 / 2 satisfies 

2 vt — 2edi v(h' £ (v 2 )\7v) + div(w) + u ■ Vu = 0, 


(2.29) 


(2.30) 


I/O 

multiplying (|2.30|) by (v-k)+ with k > ||u(-, 0)||z,°o(0) = ||po|| L oo( n) and integrating by 
parts yield 


d_ 

dt 


J (v — k) 2 + dx + 2ae J v 2a 2 |V(u — k)+\ 2 dx 


< C 


v 4 2a \u\ 2 dx 


' A k (t) 


+ ae J v 2a ~ 2 \V(v - k)+\ 2 dx, 


(2.31) 


where A k (t) = {x € I2|u(x,t) > k}. It thus follows from (|2.29D and Holder’s inequality 
that 


v 4 2a \u\ 2 dx 


JA k (t) 

,c:( 

sc( 

<c( 


v 2 \u 


‘A k {t) 


\ 2 /( 2 +«) 
1 2+£ dx J 


e/(2+e) 


^ (4+4e—2(2 +e)a)/e 




' A k (t) 


’A k (t) 


(/9 4 ( a+1 ) E 1 4 - p- 4 (a+i)e 1 


e/(2+e) 


)dx 


(2.32) 


\dA k (t) 

<C\A k (t)\ £ /\ 


(p £ ~ 2 + p~ £ ~ 2 )dx 


£(4-e)/(6(2+e)) 


|4fc(t)| 


e/6 


where (|2.27|) has been used in the last inequality. Putting (I2.32P into (12.311) leads to 


I' k {t)+ae J p a 4 \V(v-k) 


— k) + \ z dx < Cu £ J & , 


(2.33) 


where 


h{t) = / (v - k)+(x,t)dx, v k = sup \A k (t)\. 
J 0 <t<T 


Since I k ( 0) = 0, without loss of generality, we can assume that there exists some a > 0 
such that 

I k {a) = sup I k {t). 
o <t<T 
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It follows from (12.331) that 


Ikip) + j P a 1 |V(v - k) + \ 2 (x,a)dx < Cu^ 6 , 
which, together with Holder’s inequality and (|2.27|) . gives 

4W + || V (V - *)+(■, a) ||| 24/(12+ e )(n) 

< Cvf + J p a ~ l \V{v- k) + \ 2 (x,a)dx (^J p l2 ^~ a ^ e {x,a)d^j ' (2.34) 

< Cvf. 

Then, for any h> k> ||i;(-,0)Hx, 00 (O), direct computations yield 

\A h (t)\(h - k ) 2 

< \\(v - k)+(-,t)\\ 2 L 2 m 

< \\(.v - k) + (-,<j)\\ 2 L2{a) 

<C||(u-fe) + (-,a)||2 24/e(Q) |H fc ( C T)| 1 - £ / 12 

< C (||(u - *)+(•, a)||| 2(n) + ||V (v - k) + (-,a) ||^ 24 /( 12 +e )(n) ) vl ~ e/12 

< CV fc 1+e/12 , 

where in the last inequality one has used (|2.34l) . This implies 

u h < C(h - k)~ 2 v l k +£/1 \ 

which, together with the De Giorgi-type lemma 001 Lemma 4.1.1], thus shows 

sup ||p|| L oo (n) < C. (2.35) 

0 <t<T 

Finally, since w = v~ l satisfies 

2wt + 2u ■ Vw — wdivu + 4:sh' e (p)w~ 1 \Vw\ 2 = 2edi v(h' £ (p)Vw), (2.36) 

multiplying (12.361) by (w — k) + with k > ||u;(-,0)||^oo= ||/9 0 l ^ 2 \\L°°(n) yields that 


dtj 

(w — k) 2 + dx + 2ea J p a 

1 \V(w — k) + \ 2 dx 

< c 

/ w\u\\Vw\dx + C 

/ (w — k) + \u\\Vw\dx 


^A-k ( t ) J 

Ak(t) 

< c 

/ p~ a \u\ 2 dx + ea / 

p a ~ x \y{w — k) + \ 2 dx, 


^ A k(t) J 



where Ak(t) = {x € U\w(x,t) > k}. It follows from Holder’s inequality, (12.29 1) . and 
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Q2jZ7D that 


’A k (t) 


p a \u\ 2 dx 


' A k(t) 


P 


-a- 2/(2+e) 


(p 1 /( 2+£ )M) 2 da; 


< C 


< c 


\ 2/(2+e) 

M*(t) / 


f [ /9 -( 2 +( 2 + £ )“) £ 1 dx) 


e/(2+e) 


(2.38) 


'hfc(t) 


p— 6(2+(2+e)a:)/(e(4—e)) ^ 


(4—e)e/(12+6e) 


l^feWI 


e/6 


< c^ /8 , 


where = sup |Hfc(t)|. Hence, putting (|2.38l) into (12.371) leads to 

0<t<T 

— J (w — k) 2 + dx + ea J p a ^ 1 \V(w — k) + \ 2 dx < CD £ J 6 . (2.39) 

Using (12.391) and (12.271) . one can proceed in the same way as the proof of (12.351) to 
obtain that there exists some positive constant C > C such that 

sup p~ 1 (x,t)<C , 

(x,t)££lx (0,T) 

which, combined with (I2.35[) . gives (12.251) and finishes the proof of Lemma 12.41 

We still need the following lemma concerning the higher order estimates on (p, u ) 
which are necessary to obtain the global strong solution to the problem (11.151) - (11.171) 


Lemma 2.5 For any p > 2, there exists some constant C depending on e,p, and T 
such that 

r T , 

IIP _L \\(„ X7 r. „,M|P 


1 0 


Pt,V pt,u t \\ p LP{n) + \\(p, Vp,«)H^ 2iP(n) J dt < C. 


(2.40) 


Proof. First, it follows from (|2.25|) . (|2.29l) . (|2.4D . and (12.91) that 

rT 


su P m (ll n lli 2 + e (st) + ll v Plli 2 (a)) + / (ll tt lli^+2e(n) + 


0 <t<T 


dt < C, (2.41) 


which, together with the standard Holder estimates for (12.301) . yields that there exist 
positive constants C and cr G (0,1) such that 


IMIc' T > CT / 2 (f2x[0,T]) — C- 

Next, it follows from (12.301) that v = p 1 / 2 satisfies 

2 vt — 2ediv(/i(.(p)Vu) = — div(m; + X7w) — / u • Vvdx, 

where for t > 0, w(-,t) is the unique solution to the following problem 

{ Aw = u ■ Vu — |f2| _1 f u ■ Wvdx, x € H, 
f wdx = 0. 


(2.42) 


(2.43) 


(2.44) 
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Applying standard L p -estimates to (|2.44l) yields that Vtu satisfies for any p > 2 


l|Vra|| L p ( Q) < C'(p)||«|| i 2 + e( n ) ||Vw|| i , ( n ) < C(p) ||Vu|| L « (n) , 


where 


- 4 - + - - — > I. 

q p 2 2 + e p 


(2.45) 

(2.46) 


Since (12.411) and (I2.25|) imply 


u ■ \7vdx 


< C||n|| L2(Q) ||Vp|| L 2 (n) < c, 


applying standard parabolic L p -estimates to (12.431) yields that for any p > 2 
J o \\Vv\\ p LP{n) dt<C(p) + C(p) {\\u\\ p LPm + \\Vw\\ p LP{ ^dt 


< 


C{p) + C(p) [ 1 

Jo 


1 


\u\ 


lp(Q) 


dt+ z II v ^IIlp(o)^> 


where in the second inequality, (12.451) . (12.461) . and ()2.41D have been used. Thus, 


£\\Vv\\ p LP(n) dt<C(p) + C(p)f^ 


I p 

I Lp(n) 1 


(2.47) 


which, together with (12.411) . gives 


£ UVull£t?f e(n) dt < C. (2.48) 

Next, note that (| 1.15 1) o implies that it satisfies 

u t - + y/e)p~ 1 h £ (p)Au - []^p~ l h e {p) + (1 + y/£)p~ l g e (p) \ Vdivu = F, (2.49) 


where 


F = —u • Vu + (^ + Vz)p 1 V/i e (p) • Vu+ ^p l Vu-Vh £ (p £ ) 

+ (l + y/e)p~ 1 V(g £ (p))divu-p^VP-e^ 3 (p~ 1+£ 2 + p~ l ~ £ 2 ) 


(2.50) 


u. 


Since \fp{= v) satisfies (12.421) . applying the standard F p -estimates to (|2.49|) (|2.50l) (12.31) 
with periodic data, we obtain after using (12.481) and (|2.4ip that 


l 2 + £ (Qx(o,t)) + l|V 2 M||L2+e(nx(o,r)) < C + C||Vu||.£,4+2e(fi X (o,T))- 


(2.51) 


It thus follows from the Sobolev inequality ( [29J Chapter II (3.15)]) that for any 
q > 0 there exists some constant C{rj) such that 

||Vlt|| L 4+2 £ ( f2x ( 0jT )) <77(||tq||L 2 + e (f2x(0,T)) + ||V"u|| L 2+ £ ( Ox ( 0i T))) 

+ C(q)\\u\\ L 2+e(fix(o, T )), 

which, together with (12.411) and (I2.5ip . gives 


L 2 +qax(o,T)) 


+ l|V 2 d 


L 2 + £ (Qx(0,r)) ^ 


< c. 
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This, combined with the Sobolev inequality ( [221 Chapter II (3.15)]), leads to 

sup |M| L oom) < C , (2.52) 

0 <t<T 

which, along with (I2.47p . shows that for any p > 2 

£ \\Vv\\ P LP{n) dt<C(p). (2.53) 

This, together with (12.5211 . (12.251) . and the standard L p -estimates of the parabolic sys¬ 
tem (I2.49P (12.501) (12.31) . yields that for any p > 4, 

ll^t ||i/p(Qx(o,T)) + l|V 2 M||i,p(nx(o,r)) 

— C(p) + C(p) || Vtl || L 2(qx(0,T))II Vu IIl°°(Qx(0,T))II ^11 L 2 P (Qx (0,T)) 

< C(p) + - ||«t||LP(nx(0,T)) + 2 ll^ 72 ' u lkp(nx(0,T))) 

where in the second inequality one has used (12.411) and the Sobolev inequality ( [ 251 
Chapter II (3.15)]). Thus, it holds that for any p > 2, 

IWI.LP(fix(o,T)) + l|V 2 ,u llLp(nx(o,r)) < C(jp). (2-54) 

With ([2.541) and (12.531) at hand, one can deduce easily from (12.301) . (12.31) . and (12.421) 
that for any p > 2 , 

IIPtllLJhO.T.VKMn)) + l|V 2 /3||Lp(o ) T,W 1 ’P(f2)) < C(p), 

which, together with (12.541) and (12.411) . gives the desired estimate (12.401) and finishes 
the proof of Lemma 12.51 

2.2 Compactness results 

Throughout this subsection, it will be always assumed that a and 7 satisfy the condi¬ 
tions listed in Theorem 0 

We first construct the initial data. Set 

no = (8(a + 7 + 2))- 8 . (2.55) 


Choose 

satisfying 


p 0£ ec°°(n), 0 < po e < e~ 4a ° 


~ Pollii(Q) + \\P0e ~ PoIIlt(Q) + ||V(poe — Pq l ^)\\L?(Cl) < £ - 


For v >2 suitably large such that v{a — 1/2) > 5, define 

,2/(i/(2t*-l)) 


p 0£ = (p^- 1/2) +e 4 -°K«- 1/2) )' 


It is easy to check that 


li™. llPoe — PolUifn) = 0 

e->0 ' 


(2.56) 

(2.57) 
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and that there exists some constant C independent of e such that m holds. 
Since po,mo satisfy () 1.8 1) . we choose wq £ € C°°(Q) such that 


{woe - ^o/Po 1+r? ° )/(2+?70) |lL 2 +™(f2) 


< £. 


Set 


UOe = 


-1/(2+970) 

POe 


W Oe- 


Then, we have 

lim \\poeUoe - m 0 || L i(n) = 0. 


(2.58) 

(2.59) 


Moreover, there exists some positive constant C independent of e such that (12.21) holds. 

Extend then ( pq £ ,uq £ ) fi-periodically to K 2 . The standard parabolic theory [29] , 
together with Lemmas 12.41 and 12.51 thus yields that the problem (11.151) (11.171) (12.31) , 
where the initial data (po,uo) is replaced by ( po £ ,uo £ ), has a unique strong solution 
(Pe,u E ) satisfying 


Pe, U £ , p et , U £t , V 2 p e , V 2 u £ € L p {yi X (0,T)), 

for any T > 0 and any p > 2. Moreover, all estimates obtained by Lemmas 12.1112.31 still 
hold for ( p £ ,u £ ). 

Letting e —>• 0 + , we will modify the compactness results due to [35] to prove that the 
limit (in some sense) {p iy fpu) of (p £ ,y/p^u £ ) (up to a subsequence) is a weak solution 
to (11.11) (| 1.3 D (11.61) (11.71) . We begin with the following strong convergence of p £ . 

Lemma 2.6 There exists a function p € L°°(0, T; L 1 (fl) n L 7 (fl)) such that up to a 
subsequence, 

p £ —>• p in L 7 (fl x (0, T)). (2.60) 

In particular, 

p £ —>• p almost everywhere in Q x (0 ,T). (2-61) 


Proof. First, it follows from (12.41) and (12.91) that there exists some generic positive 
constant C independent of e and T such that 


sup / (p e \u e 2 + p £ + pf)dx + 

[ [ h £ (p £ )\Vu £ 

0 <t<T J 

Jo J 

+ £ Jo J Pe + 

u £ \ 2 )dxdt 


: i £ \ 2 dxdt < C, 


2 dxdt 


(2.62) 


and that 


sup [ p £ l (ti £ (p £ )) 2 \S7p £ \ 2 dx + 
0 <t<T J 



pf^-^V p £ \ 2 dxdt 


+ £ 


13/3g—£“ 


sup 

0 <t<T 


f £ 2+7 " 1 + pf e 2 " 1 / 8 


dx < C. 


(2.63) 


Then, (|2.62l) and (|2.63D imply that 


sup ||Vp"|| i27 /( 7 +i)(o) < C sup ||p 2/2 || L 2 7(n) sup ||Vp“ 1/2 || L*(p.)<C. (2.64) 

0 <t<T 0 <t<T 0 <t<T 
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Moreover, note that p" satisfies 

(p“)t + div(p £ u £ ) + (a - l)p"divu £ 

= eadiv(p“ _1 /i'(p £ )Vp £ ) - ea(a - ^)p £ _2 d £ (p £ )|Vp £ | 2 


It follows from (I2.62|) . (12.631) . and (I2.24P that 

sup ||p“u E || L i(n) < C sup ||p"^ 1/2 || L 2 (n) sup \\pl /2 u £ \\ L 2 { n) <C, 


0 <t<T 
r T 


0 <t<T 


0 <t<T 


llp"divu £ || 2 1(n) df < C ||p £ /2 || 2 2 ( q) ||p £ /2 divu E || 2 2 (f2) dt < C, 


and that 


sup / (p“ 1 h' £ (p e )\Vp e \ + p £ 2 /i^(p £ )|Vp £ | 2 ) dx 

' ;<T./ 

c sup f (p 2a_1 + p7 1 (^(p £ )) 2 |Vp £ | 2 ) dx < C. 
0 <t<T J 


0 <t<T 
< 


The combination of (12.65|) — (I2.68[) implies that 


(2.65) 


( 2 . 66 ) 

(2.67) 


( 2 . 68 ) 


II (Pe )t llL 2 (0,T;W /_1 ’ 1 (f2)) < C. (2.69) 

Letting e —>• 0 + , it follows from (12.641) . (12.691) . and the Aubin-Lions lemma that up to 
a subsequence 

p“^p“ inC([0,T];L 3 / 2 (L!)), 


which implies that (12.611) holds. In particular, it holds that 


pf~ 1 / 2 p"" 1 / 2 in L 2 (0,T;L 2 (L?)). 


(2.70) 


Finally, it follows from the Sobolev inequality, (12.621) . and (|2.63l) that 

[ T Wn'y ||( 5 7+3(a-l))/(37) r U 
/ HPe ll L (5 7 +3(a-l))/(3 

•J 0 

< c £ IklSln, (llftllHS 1 + l|v 4 7 + “- 1 )/ 2 lly (tl) ) it < C, 


(2.71) 


which together with (12.611) thus gives (12.601) due to (57 + 3(a— l))/(3y) > 1. The proof 
of Lemma 12.61 is finished. 

Before proving the strong convergence of \J~pfu e in L 2 (fl x (0, T)), we show first the 
following compactness of p £ 7+1 ^ 2- u £ . 


Lemma 2.7 There exists a function m(x,t) G L 2 (fl x (0,T)) such that up to a subse¬ 
quence, 

p( 7+1)/ V -+min L 2 { 0, T; L p (0)), (2.72) 

for all p G [1,2). Moreover, 

p(7+!)/ 2 U£ p('y+ 1 i/ 2 u a l mos t everywhere (x,t) € Tl x (0,T), (2.73) 


where 


u(x,t ) = 


m(x, t)/p( 7+1 )/ 2 (x,t) 

0 , 


for p(x, t ) > 0, 
for p(x, t) = 0. 


(2.74) 
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Proof. First, since a € (1/2, (1 + 7)/2], it follows from (12.621) . (12.631) . and (12.241) that 
for any rj > 0, 


j o l|V(4 1+7)/2 M £ )||ii (n) dt 
< cj_' 11 /3 ( 1 +'r-a )/2 111 2 (n) 11 p“ /2 Vu £ 111 2 (n) dt 
+G^ T |i4 7 ‘ 2Q+2)/2 ^iii 2( f 2) iivpr 1/2 ii| 2(f2 )dt 

y y 

<C{rj) + C(rj) \\pI /2 u £ \\ 

L 2 (fi)^ + C V J o \\p£ +l)/2 Ue\\ 
which together with the Sobolev inequality gives 

J q \\V(pi 1+l)/2 u £ )\\ 2 Ll(n) dt +\\p ( £ 1+ ^ /2 u £ \\ 2 L2{n) dt<C, 
due to the following simple fact 

sup f p^/ 2 \u e \dx<C sup llp/ll^m sup \\pl /2 u £ \\ L 2 {n) < C. 
0 <t<T J 0<t<T K ' 0 <t<T 


(2.75) 


(2.76) 


(2.77) 


Next, we claim that 

ll(pi 7+1 )/ 2 ^E)t||Li(o,T;W—M(fi)) < C, (2.78) 

which, combined with (|2.76l) and the Aubin-Lions lemma, yields that there exists a 
function m(x,t) € L 2 (Pl X (0,T)) such that up to a subsequence, (12.721) holds for all 
p € [1,2). In particular, 

p(7+ 1 )/ 2 Ue —> m almost everywhere (x,£) gllx (0,T). (2.79) 

Moreover, since pl^ 2 u e is bounded in L°°(0, T;+ 2 (fl)), Fatou’s lemma gives 



lim inf 


lo (7+1)/2 u I 2 

— - 7y—-—dxdt < oo, 

Pi 


which implies m(x,t ) = 0 almost everywhere in {(x,t) € x (0, T)|p(x, £) = 0}. Hence, 
for u(x,t) as in f|2.74[) . we arrive at 

m(x,t ) = p^ +1 ^ 2 (x,t)u(x,t), 


which together with (|2.79p gives (12.731) . 

Finally, it remains to prove (1 2.78 j) . In fact, note that 

(4 7+1)/ V) t = ^-i / 9/ 7 " 1)/2 ( / 9 £ ) t rt £ + 4 7+1)/2 (u £ ) t . 

One can use (j 1.15 D to get 

P ( e' y ~ 1)/2 (ps)tu £ = -p < £ / - 1)/2 div(p £ u £ )u £ +£p] /2 div(pf 1/2 ti £ {p e )\7p £ )u £ 

= -^—div(4 7+1)/2 u £ (8) u £ ) - ip/ 7+1)/2 divu £ u £ 

7+1 7+1 

4-^—p( 7+1)/2 tt £ • Vu £ + ediv(4 7_1)/2 /i(.(p £ )V / o £ ® u £ ) 

7 + 1 

- y/4 7 - 3)/2 /4(p e )|Vp e | 2 u e - £ 4 7 - 1 )/ 2 / l /(p £ )Vp £ • Vu e , 


(2.80) 


(2.81) 
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and 


< 


< 


(2.85) 


p{ 'y+ 1 ')/ 2 (u e ) t = - 4 7+1)/2 «e • Vu £ + div(4 7 l)/2 h £ (p e )(Vu + y/eS/Ue)) 

- ^p^- 3)/2 h E (p E )Vp e ■ {Vu + y/eVu £ )) 

+ (1 + v^)V(pi 7_1)/2 c/ £ (p £ )divu £ ) (2.82) 

- —— ^ + ^ pi 7 ~ 3)/2 g £ (p £ )Vp £ divu £ 

- 4 7_1)/2 Vp 7 - e- £ ~\pf' 2 +( 7 “ 1 )/ 2 + p- £ ' 2 +( 7 - 1 )/ 2 ) W£ . 

One needs to estimate each term on the righthand side of (12.811) and (12.821) . It follows 
from the Holder inequality, (12.621) . (12.631) . (12.761) . and the Sobolev inequality that 

[ T \\p^/ 2 \u £ \ 2 \\ LHn) dt < C [ T II(p £ + P 2 +1 )\u £ \ 2 \\ LHn) dt < C, (2.83) 

Jo Jo 

[ ||4 7+1)/2 K||V« £ ||| L i (n) cft 

0 T T (2.84) 

<cJ q ii(^ 2 + 4 7+1)/2 )klll i 2 (n) 

£ / ||pi 7_1)/2 /i £ (p £ )|Vp £ ||u £ ||| L i ( n)dt 

Jo 

Ce j J p7 1 /i £ (p e )|Vp £ | 2 |rt £ | 2 dxdt + C j J p 7 ' 1 h £ (p £ )dxdt 

C + C Jo I (P " +a_1 + + p ?~ 5/ ^ dt ~ C ’ 

£ [ IIPe 7_3)/2 ^(Pe)|Vp £ | 2 |rt £ ||| Z/ i ( n)dt 
Jo 

<CeJ o j (p- 1 h'(p £ )|Vp £ | 2 |u £ | 2 + p 7 - 3 h £ (p £ )|Vp £ | 2 )dxdt<C, 

J p^ 1]/2 (h £ (p £ ) + \g £ (p £ )\)\Vu £ \dxdt 
< C J J h £ (p £ )\Vu £ \ 2 dxdt + C J J p]~ 1 h £ (p £ )dxdt < C, 

J o j (pi y ~ 1)/2 K(pe) + Pp~ S)/2 ( h e(p £ ) + \9 E (p E )\))\Vp £ \\Vu £ \dxdt 

<cjj J p[ 7_3)/2 h £ (p £ ) | Vp e 11 Vti £ | dxdt 
< C j J h £ (p £ )\S7 u £ \ 2 dxdt + C J J p]~ 3 h e (p e )\V p £ \ 2 dxdt < C, 

f J Pe 7_1)/2 |Vp 7 Mxdt 

<c [ 7 \\pr a/2 \\LW\Vp^ +a - 1)/2 \\ L * m dt<C, 

Jo 


( 2 . 86 ) 


(2.87) 


( 2 . 88 ) 


(2.89) 
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< Ce 


+ Ce 


£ j (pf 2+(7 " 1)/2 + p~ £ ~ 2+{l - l)/2 )\u £ \dxdt 

f J{P £ s 2 +P7 E 2 )\u e \ 2 dxdt 
1 £ J (pf 2 ^- 1 + P r~ 2+i - l )dxdt 


<C + Ce 


i o 


(P\ 


f 2 +7-1 i „-e 


+ £ )cfa;dt < (7, 


where in the last inequality one has used the following simple facts that 

rT 


£ 3 J j p e e 2+7 1 dxdt 

e -£ _1 /(6(l+(7-l)e 2 )) 

e 13(l+( 7 -l)£ 2 )/(3(l+(7-l) £ 2 )) 


(7e- 1/(9£ ) / 13/3 _- 3 f 

< -- sup e 13 ' 3 e / 

£ 3 0<t<T \ J 

< C'e _13 e _1/(9 ^ -+ 0 as e -+ 0, 


(e 


13/3 -e 3 e 2 +7-l'|f=7 


) e 2 +7—1 dxdt 


e~ +"/—l 

Pl^-'dx' 7 ^ 


and that 


p £ dxdt 


e -s 1 /(8+e 2 ) j'T 
e 104/(3(8+ £ 2)) J Q 
C' e - 1 /(9e) 


£r 13 / 3 e _e_3 p _£ - 2_1/8 ^ 




< —^- sup e 13 / 3 e _£_3 / 

£ 5 o<t<T V J 

< Ce~ b e~ l ^ 9e ^ —> 0 as e -+ 0. 


- 1 / 8 * 


8/(8+e 2 ) 

dxdt 
8/(8+e 2 ) 


(2.90) 


(2.91) 


(2.92) 


Thus, all these estimates (I2.83p - (|2.9UD together with (I2.80p - (I2.82D yield (12.781) . The 
proof of Lemma 12.71 is completed. 

Now we are in a position to prove the strong convergence of yfpfu e in L 2 (Q x (0, T)) 
which in fact is essential to obtain the existence of global weak solution to the problem 

CCD (USD CLiID (H2D- 


Lemma 2.8 Up to a subsequence, 

yfpfu £ —> yfpu strongly in L 2 (0, T; L 2 (fl)), (2.93) 

with 

y/pu € L°°(0, T; L 2 (0)). (2.94) 

Proof. First, Lemma 12.31 yields that there exists some constant C independent of e 
such that 

sup / p £ \u £ \ 2 ln(e + \u £ \ 2 )dx < C, (2.95) 

o <t<TJ 
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which, together with (12.730 . (12.610 . and Fatou’s lemma, gives 

J j p\u\ 2 ln(e + \u\ 2 )dxdt < a liminf p £ \u £ \ 2 ln(e + \u e \ 2 )dxdt 


< 


liminf [ [ p £ \u e \ 2 ln(e + \u e \ 2 )dxdt < C. 

£ ^° Jo J 


Next, direct calculation shows that for any Ad > 0, 

rT 


(2.96) 


II | 'fpju £ - y/ i pu\ 2 dxdt 

< 2 / / \Vfeu e lQ Ue \<M) - VP ul (\u\<M)\ 2 dxdt (2.97) 

Hi \y/~Pe u e^-{\ue\>M) I dxdt + 2 II Wpul { \u\>M)?dxdt. 


Next, it follows from (12.730 and (12.610 that yfpju e converges almost everywhere to 
yfpu in the set {(x,t) E J) x (0,T)\p(x,t) > 0}. Moreover, since 


y/Jh\ u e\^-{\u e \<M) — ^-\fPei (2.98) 

and p £ —>• 0 almost everywhere in the set {(x,t) € O x (0,T)\p(x,t) = 0}, we have 
yfpl u e\\u e \<M) VP ul (\u\<M) almost everywhere in O x (0,T), 
which, together with (12.98(1 and (12.600 . implies 

II \\fPe u e^-(\u e \<M) - VP u ^-{\u\<M) I dxdt ->• 0. 

Next, it follows from (12.950 and (12.960 that 


\\fp-£ u e l {\ue\>M)\ 2 dxdt + 


Wpu\\ u \>M)\ 2 dxdt 


< 


< 


ln(e + 
C 


ln(e + M 2 ) 

Substituting (12.990 and (12.1000 into (12.970 yields that up to a subsequence 

,T , c 


limsup / 

£—^0 J 0 


| y/lhU e - pu\ dxdt < 


ln(e + M 2 ) 


(2.99) 


f M 2 ) J j (pe\ u e? + \u e \ 2 ) + p\u\ 2 \n(e + \u\ 2 )) dxdt (2.100) 


( 2 . 101 ) 


for any M > 0. We thus obtain (|2.93Q by taking M —>• oo in (12.1010 . 

Finally, the combination of (12.620 with (12.930 gives (|2.940 immediately. The proof 
of Lemma 12.81 is completed. 

As a consequence of Lemmas 12.61 and 12.81 the following convergence of the diffusion 
terms holds. 
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Lemma 2.9 Up to a subsequence, 


p £ Vu £ -> p°X7u in V , (2.102) 

p“(Vu £ ) tr -> p a {Vu) tr in V, (2.103) 

p“divu £ -> p“divw m P'. (2.104) 

Proof. Let ^ be a test function. Then it follows from (12.631) . (12.701) . and (|2.93l) that 

J j p £ X7u £ (pdxdt 

= - j J p ( f~ 1/2 y/p~ £ u £ S7(j)dxdt - J J V Pe~ 1/2 \fpl u e<t><lxdt 

r T p c\ C T p 

—i — J J p a ~ 1 ^ 2 yPpu\/(j)dxdt — -— j- J j Vp“^ 1//2 yfpu<j>dxdt, 


which gives (|2.102l) . Similar arguments prove (12.1031) and (12.104p . and the proof of 
Lemma 12.91 is completed. 

2.3 Proof of Theorem ll.lt = T 2 

First, rewrite (|1.15|) i as 


{p £ )t + di v(p £ u £ ) = ediv(h' £ (p £ )Vp £ ) - -p £ 1 h' e (p £ )\Vp e \ 


(2.105) 


It follows from (I2.62|) and (12.631) that 

rT 


J J K(p £ )\Vp £ \dxdt 


f f T f \ 1 / 2 

< Ce yj J P7 (K(p E )) 2 \Vp £ \ 2 dxdt 

< Ce, 


(2.106) 


and that 


J J Pe 1 K(Pe)\^ / Pe 1 2 dxdt 

fCe(^J jpf l (ti e (p £ )) 2 \Vp £ \ 2 dxdt^j (^j Jp- 1 \Vp e \ 2 dxdt S J 


1/2 or 

< Ce 2 / 3 (e 2 / 3 £ j (p- 5 / 4 + p 2 p~ 8 / 3 ) \Vp £ \ 2 dxdt ) 

< Ce 2 ' 3 . 


1/2 


(2.107) 


1/2 


Then, letting ip be a test function, multiplying (|2.105D by ip, integrating the resulting 
equality over P x (0, T), and taking e —>• 0 (up to a subsequence), one can verify easily 


after using (12.601) . (12.931) . (12.571) . (12.1061) . and (12.1071) that ( p,yfpu ) satisfies (11.101) . 
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Next, (11.151) implies that 

(p e u £ )t + di v{p £ u £ <g> u £ ) - di v(p/Vu £ ) - (a - l)V(p"dmi £ ) + VP(p £ ) 

= edi v(ti e (p £ )'Vp £ <g> u £ ) - ^p^ti^p^lVp £ \ 2 u £ - eh' £ (p £ )V p £ ■ Vu £ 

-e~ £ \p e £ 2 + p~ £ 2 )u £ + y/ediv(h £ (p £ )Vu £ ) + y/eV (g £ (p £ )divu £ ) 

e-1/3 

+ e 1 / 3 di v((pj/ 8 + pl)Vu £ ) -— V(/9^ 8 div« e ) + £ l ^(j — l)V(pJdivu e ). 


(2.108) 


Using (|2.62l) and (|2.63l) . we have 

£ l / h e(Ps)\^Pe\\ u e\ dxdt 

<Cs(^J j pj 1 (h! e (p e )) 2 \V p e \ 2 dxdt) ( [ [ p e \u e \ 2 dxdt) 


J \Jo 


<Ce^ 0 , 




(2.109) 


Pe 1 ti E ( yPe)\V Pe\ 2 \u e \dxdt 



1/2 


P £ h £ (p £ )\\7p £ \ \u £ \ dxdtj 
1/2 

P^KiPeW Pe?dxdt 


where (|2.107l) has been used in the last inequality, 
h' £ (p £ )\S7 p £ \\S7 u £ \dxdt 



< Cy/e 

< Cy/e-> 0, 


( 2 . 110 ) 


(e/ J P £ 1 h' £ {p £ )\\/p £ \ 2 dxdt S j (^j J h £ (p E )\Vu £ \ 2 dxdt^ 


( 2 . 111 ) 


/o 


(Pe 2 +P £ e 2 )\u £ \dxdt 


_ r T r _ _ \ 1 / 2 

< C ( e _£ 3 J J (p £ e 2 + p“ £ 2 )|tt e | 2 dxdt 


x e 


< C e 




(Ps 2 + Pe £ 2 )dxdt j 
'o J J 


1/2 



(Pe + Pe 2+7 1 +Pe £ 


1/2 


due to (12.91 [) and (|2.92l) . 


( 2 . 112 ) 


VSJ J (h £ (p £ ) + \g £ (p £ )\)\Vu £ \dxdt 

/ rT r \ 1/2 / rT r \ 1/2 

< Cy/e ( / / h £ (p £ )dxdtj ( / /i £ (p e )|Vtt £ | 2 dxdt 

< Cx/e ->■ 0, 


(2.113) 
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where in the second inequality one has used the fact that 


J o H/ 9 ellLT(n) c ^ ~ J o H^ e lll7+i/6(n)^ 

J 0 (HaIIl 1 ^) + 1 1 y°e 1 1 ^5^/s+l- 1 (^2) ) dt 

<C, 


due to 7 > 1, a > 1/2, and (12.7111 . 




h e (Pe) |Vu e 



1/2 


(2.114) 


(2.115) 


where (I2.114|) has been used in the second inequality. 

Finally, let 0 be a test function. Multiplying (12.10811 by </>, integrating the resulting 
equality over 11 x (0,T), and taking e —» 0 (up to a subsequence), by Lemmas 12.61 
12.81 and 12.91 we obtain after using (12. 1091) — (12 .1131) . ()2. 1 15[) . and (12.5911 that {p,yfpu) 
satisfies (| 1.11 [) . The proof of Theorem 11.11 in the case 9 = T 2 is completed. 


3 Proof of Theorem 11.3b Vt = T 3 

In this section, we will show how to modify the analysis in the previous section to deal 
with the 3-dinrensional case with periodic boundary conditions. 

3.1 A priori estimates 

For e as in (11.1711 . let 9 = T 3 and smooth functions po £ > 0 and uq £ satisfy (12.111 and 
(12.211 for some constant C independent of e. Moreover, if a € (1,2), in addition to (12.lH 
and (12. 2 j) . it holds that for some C independent of e, 

J P0e\u0e\ 4 dx < C. (3.1) 

We extend po e and uo £ 9-periodically to M 3 . For 7 as in (11.1711 . consider 

Pt + div (pu) = e p 1 / 2 div (p- 1 h' e (p)\7 p), 

< put + pu ■ Vu — div(h e (p)Vu) — \7(g £ (p)divu) + VP (3-2) 

=-e~ e ~\p e ~ 2 + p~ £ ~ 2 )u, 

with h e (p) and g £ (p) as in (11.1611 . The initial condition for the system (|3.2j) is imposed 
as: 

(p,u)(x, 0) = (po £l u 0£ ), x € 9. (3.3) 

Let T > 0 be a fixed time and (p, u ) be the smooth solution to (13.211 (11.1611 (11.1711 (13. 3 p 
on 9 x (0, T\. 

After some minor modifications, one can check easily that all the estimates in Lem¬ 
mas [2J] and [272] still hold for cc € [3/4,2) and 7 € (1,3). That is 
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Lemma 3.1 Let a E [3/4,2) and 7 € (1,3). Then there exists some generic constant 
C independent of e and T such that m and (12.91) hold with Ll = T 3 . 

To obtain the Mellet-Vasseur type estimates for the three-dimensional case, we need 
to impose some additional constraints on 7 and a. 

Lemma 3.2 Assume that a E [3/4, 2) and that 7 E (1,3) satisfies 7 E ((a + l)/2, 6 a — 
3) in addition. Then there exists some generic constant C depending on T but inde¬ 
pendent of s such that 


sup / p(e + |u| 2 ) ln(e + \u\ 2 )dx < C. (3.4) 

0 <t<T J 


Proof. It is easy to check that (12.221) still holds. Hence, it remains to estimate the 
righthand side of (12.221) . In fact, since 7 E ((a + l)/2, 6a — 3), we have 6 = € 

(0,1). Then 


Jo 

< 

< 


cf (Jpln 2/ ‘( e + \u\ 2 )dx\ (^J dt 

C + C lo I pa ~ 1+5l/3dxdt - C ’ 


(3.5) 


where in the last inequality one has used (12.711) . This together with (12.221) yields (13.4|) . 
The proof of Lemma 13.21 is completed. 

Furthermore, the following estimates on the L°°(0, T; L 4 (H))-norm of p x ^u will be 
used later. 


Lemma 3.3 Assume that a E [3/4,2) and that 7 E (1,3). Then there exists some 
constant C(e) depending on e and T such that 

p\u\ i dx + J J h e (p)\u\ 2 \Vu\ 2 dxdt < C(e). (3-6) 

Moreover, if in addition (a + l )/2 < 7 < 3a — 1, there exists some constant C\ inde¬ 
pendent of s such that 

p\u\ 4 dx + a h £ (p)\u\ 2 \SJu\ 2 dxdt < C\. (3.7) 


sup 
0 <t<T. 


sup 
0 <t<T. 
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Proof. First, multiplying (13.21) 9 by \u\ 2 u and integrating by parts give 
p\u\ 4 dx)t + j h e (p)\u\ 2 \Vu\ 2 dx 4 - 2 J h e (p)\u\ 2 \V\u\\ 2 dx 
+ / g £ (p)(divu) 2 \u\ 2 dx + 2 / g £ (p)divu\u\u ■ V\u\dx 


1 

4 ( 


+ e 

e 

4 


J(p £ 2 + p~ e 2 )\u\ 4 dx 

■ j p 1 / 2 dw(p~ 1 / 2 h , £ (p)'Vp)\u\ 4 dx + J Pdiv(\u\ 2 u)dx 
= — ^ J p- l h' e (p)\S7p\ 2 \u\ A dx + e J h £ (p)\u\ 3 V p ■ V|u|dx 
+ J Pd\v(\u\ 2 u)dx 

< 4e j ph' £ (p)\u\ 2 \Vu\ 2 dx + J Pdiv(|u| 2 u)dx 

< 4(a + 7 )e J h £ (p)\u\ 2 \S7u\ 2 dx + J Pdiv(|u| 2 u)dx. 

Next, Cauchy’s inequality implies that for any a € [3/4,16/5], 

— 2 J p & \u\ 2 \V\u\\ 2 dx + (1 - d) J p“(divu) 2 |u| 2 da; 

+ 2(1 — a) J p“divu|u|u • V|u|dx 

< (1 — a + ^ ^ f p a (divu) 2 |u| 2 dx 


< — [ p & (dwu) 2 \u\ 2 dx 
32 / 


< 


P J p Q |Vu| 2 |u| 2 dx, 

which, combined with (11.161) . shows 

— 2 J h £ (p)\u\ 2 \X7\u\\ 2 dx — J g e (p)(divu) 2 \u\ 2 dx 
! J g £ {p)d\vu\u\u • V|u|dx 


- 2 


< 


27 

32 


J /i e ( / o)|u| 2 |Vu| 2 dx. 


(3.8) 


Substituting this into (|3.8j) yields that there exists some constant C\ independent of 
e such that 


sup / p\u\ 4 dx + / / h £ (p)\u\ 2 \\7u\ 2 dxdt 

o <t<T J Jo J 

< Ci + Ci [ T V2 

JO 


(3.9) 


dt, 
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where one has used the following estimate 


j P\div(\u\ 2 u)\dx 

p a \u\ 2 \Vu\ 2 dx + C\ J p 2l ~ a \u\ 2 dx 
4/p1n| 2 |V.| 2 dx + Ci||p 2 —/ 2 || L2(n) (l + ||p 1/4 w||^ 4(n) 



Then, if 7 € [(a + l)/2, 3a — 1], it holds that 


£ \\p 2 ^- a - 1/2 \\ L 2 m dt < C, £ (||p||^ 2 (n) + \W +2Q - 3/2 \\mn)) dt < C U (3.10) 
where in the second inequality one has used 


f 


|| p 7 + 2 Q -3/2|| L2(o)dt < 

Jo 

< Ci + Ci 


«-1/2 | 
I,6“-3(Q)I 




/ T ||Vp^+' 

J 0 


a—1)/2 112 


L 2 (Q)^ — Cl 


(3.11) 


due to (12.41) . (12.91) . and the Sobolev inequality. It follows from (13.91) . (13.101) . and the 
Gronwall inequality that (13.71) holds. 

Finally, it follows from (12.91) that 


sup ||p 27 " 1/ -\\ L 2 (n) < c(e), 

0 <t<T 


which together with (13.91) and the Gronwall inequality gives (13.61) . The proof of Lemma 
13.31 is thus completed. 

With (12.41) . (12.91) . and (13.61) at hand, similar to Lemmas 12.41 and 12.51 one can prove 


Lemma 3.4 Assume that a E [3/4,2) and that 7 € (1,3). Then there exists some 
constant C depending on e and T such that for all (x,t) € Q x (0, T) 

C- 1 < p(x,t ) < C. 

Moreover, for any p > 2 , there exists some positive constant C depending on e,p, and 
T such that 

J o (ll(P) u )t\\ P LP {ci) + ll(P.“)llw 2 ,p( n) ) dt < C. 


3.2 Compactness results 

Throughout this subsection, it will be always assumed that a and 7 satisfy the condi¬ 
tions listed in Theorem [L31 

We first construct the initial data. Choose po )£ as in ([2.561) . Hence, (|2.57[) and (12.11) 
also hold. If a € [3/4,1], define uq £ as in (|2.58[) . If a E (1, 2), set 

uoe = Po e 1 / 4 wo £ , (3.12) 
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where wq £ G C' 0O (ri) satisfies 

Ikoe - ™o/Po /4 |U*(n) < £■ 

It is easy to check that (|2.59l) and (12.21) are still valid for uq £ as in (]2.58j) . Moreover, 
(12.591) . (12.2 j) . and (13.11) hold for uo £ as in (13.121) . 

We then extend {pq £ ,uq e ) 9-periodically to R 3 . Similar to the two-dimensional case, 
using Lemma f3.4l and the standard parabolic theory [29], one can show that the problem 
m (ebd (fLiTD m has a unique strong solution ( p £ ,u £ ) satisfying for any T > 0 
and any p > 2, 


p £ , u £ , p £t , u st , V 2 p £ , X7 2 u £ G L p (Q x (0, T)). 

Lemma f3.1l thus shows that there exists some generic positive constant C independent 
of e and T such that (12.621) and (12.631) still hold with 9 = T 3 . Hence, the combination 
of (|2.62l) with (|2.63l) implies that (12.911) and (12.921) are still valid for Q = T 3 . 

Moreover, for e —>• 0 -1- , it is easy to check that Lemma [2.61 holds also for the case that 

n = t 3 . 

The following lemma deals with the compactness of the momentum. 

Lemma 3.5 If a € [3/4,1], there exists a function m(x,t) G L 2 (0, T; L 3 / 2 (9)) such 
that up to a subsequence, 

p £ u £ ^min L 2 (0,T;L p (9)), (3.13) 

for all p G [1,3/2). Moreover, 

p £ u £ —>• pu almost everywhere (x,t) G fi x (0,T), (3.14) 


where 


u(x,t ) = 


m(x,t)/p(x,t) 

0 , 


for p(x,t) > 0, 
for p(x, t ) = 0. 


Proof. First, since a G [3/4,1], it follows from (12.621) . (12.631) . and the Sobolev inequality 
that 

J \\^ {Pe u e)\\L 1 {Q)^ 

<c f (liftlii,;?», n) ilp? /2 v««iii I(n) + iiMiv Pe |i|2 1(n) ) n 

J 0 

T 

<cj 0 (\\p e \\ LH n) + l|Vpr 1/2 |lK ( ( oT 1) ) 2 ^ \\Pe /2 Vu £ \\l* {Q) dt (3.15) 

+ CJ ||kl|Vp £ ||| 2 i( n) dt 
< C + C J |||^||Vp £ ||||i (n) dL 
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Thus, if 7 > 3a — 1, then 1 € [2a — 1, a + 7 — 1] due to a € [3/4,1], Hence, it follows 
from (I2.62p and ([2.631) that 


[ lll«e||Vp e 

JO 




dt 


< 


\\p\ /2 ue\\h m (l|Vpr 1 / 2 ||! 2 ( n) + HV^ +a - 1)/2 || 2 L2(n) ) dt < C. 


(3.16) 


For 1 < 7 < 3a — 1, (12.621) . (12.63ft . (13.71) . and (13.111) imply that 

J o \\K\\VPe\\\ 2 L i ( n) dt < J o ll^ +5/4 ll! 4 (n ) II^Vlli4 ( o ) l|Vpr 1/2 |l! 2( n ) ^ 


< 


c l (HpIIlmo) + IIp 7 + 2 “ _ 3 / 2 IIl 2 ( q)) dt < c, 


(3.17) 


where in the second inequality one has used 7 + 4a > 4 due to a € [3/4,1] and 7 > 1. 
Putting ([3.16ft and (13.171) into ()3.15l) leads to 


J ||V(p e 'U £ )|| 2 i^dt — C- 

Next, similar to (12.1081) . it holds that 

[p e u £ )t + div(p £ u e <g> u e ) - div(p £ Vrt £ ) - (a - l)V(p £ divu £ ) + Vp/ 
= ediv(/i(.(p £ )Vp £ ® u e ) - ^p“ 1 /i £ (p £ )|Vp £ | 2 tt £ - sh' e (p e )V p e • \7u £ 


(3.18) 


-e £ \p e £ 2 + p e £ 2 )u e + e 1 / 3 div((pj /8 + p])Vu e ) 
,1/3 


(3.19) 


8 


V(pJ /s divu £ ) + e 1/3 (7 - l)V(p/div« £ 


Then, each term on the righthand side of (13.191) can be estimated similarly as those 
of (12.109ft - (12.1131) and (12.1151) . Moreover, for the terms on the left hand side of (13.191) . 
we have 

J J p £ |u £ | 2 dxdi + J j pjdxdt < C , 

[ [ p £ |V u e \dxdt < C f f h £ (p £ )\V u e | 2 dxdt + C f f h £ (p £ )dxdt 

Jo J Jo J Jo J 

< C. 

Hence, 

II (Pe^e)t llz, 1 ( 0 ,T;W /_ 1 ’ 1 (f 2 )) < C. (3.20) 

With ([3.18ft and (13.20[) at hand, one can finish the proof of Lemma [3.51 similarly as that 
of Lemma P 

When a € (1,2), the following compactness result of pi 1+1 ^ 2 u £ is needed. 


Lemma 3.6 Assume that a € (1,2) and that 

7 € [2a — 1, 3a — 1]. 


(3.21) 
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Then there exists a function m(x,t ) € L 2 (fi x (0, T)) such that up to a subsequence, 

pi 7+1)/2 u £ -> m m L 2 (0, T; L p (fl)), (3.22) 

for all p G [1,2). Moreover, 

p^ +l )/ 2 u e —»■ p( 7+1 )/ 2 u almost everywhere (. x,t ) € 12 x (0,T), (3.23) 


where 


u(x, t) = 


m(x,t)/ p ( 7+1 )/ 2 (x, t) 

0 , 


/or p(x,t) > 0 , 
/or p(x,t) = 0 . 


Proof. First, it follows from (I2.62p . (12.631) . and (13.71) that 
f o l|V(4 1 + 7 )/ 2 u £ )|| 2 6 / 5(n) dt 

- c Jo ll4 1+7_a)/2 Ili3 ( fi)llp“ /2v ^lli 2( Q)^ (3>24) 

+ C J 0 11 ^ /4 111 12 (O) 11 Pe /4 “e 11(O) 11'VpZ /2 11 ia (n) '* 

<<?, 

where in the second inequality one has used the Sobolev inequality, a € (1,2), and 
7 G [2a — 1,3a — 1]. Thus, the combination of (13.241) . (12.771) . and the Sobolev inequality 
shows 

£ ||V(4 1 + 7 )/ 2 u £ )|| 2 6/5(n)( it + £ ||pi 1 + 7 )/ V||| 2(n )Cft < a (3.25) 
Next, note that (12.801) and (12.811) both still hold. Moreover, it follows from (13.2 j) that 
pi 7 + 1 )/ 2 (w £ )t = - p^ +V > /2 u e • Vu £ + div(p^ 7 _ 1 )/ ' 2 h e (p £ )Vu £ ) 

- ^pi 7 - 3)/2 h £ (p £ )Vp £ • Vu e 

z _ (3.26) 

+ V(p[ 7 _ 1 )/ 2 p e (p £ )divu £ ) - p( 7 _ 3 )/ 2 p £ (p £ )Vp £ divu £ 

- p^ 1 )/ 2 Vp 7 - e- £ “ 3 ( P r 2+(7 “ 1)/2 + Pf £ ~ 2+ ^~ 1)/2 )u £ . 

Using (12.801) . (12.811) . and (13.261) . one can prove Lemma 13.61 in a similar way as that of 
Lemma 12.71 

Next, as a consequence of Lemmas 12.6113.2113.51 and 13.61 similar to Lemma [27771 one 
can obtain 


Lemma 3.7 Assume that a and 7 satisfy the conditions listed in Theorem 1 1. 31 Then 
up to a subsequence, 

yfpfu £ —> yfpu strongly in L 2 (Q x (0, T)), (3.27) 


with 


y/pu£ L°°(0,T;L 2 (Q)). 


(3.28) 
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Finally, similar to Lemma 12.91 one can prove 

Lemma 3.8 Assume that a and 7 satisfy the conditions listed in Theorem \1.3l Then 
up to a subsequence, 

p E Vu e —> p a Vu in V ', (3.29) 

p“(Vu £ ) tr ->• p“(Vu) tr in V, (3.30) 

p“divTt £ —>■ p a divu in V'. (3.31) 

3.3 Proof of Theorem 11.31 £1 = T 3 

Using (|2.62l) . (12.631) . (12.911) . (I2.92|) . and Lemmas 12.6113.71 and 13.81 one can finish the 
proof of Theorem 11.31 where £l = T 3 in a similar way as that for the 2-dimensional 
periodic case in Theorem ll.il 


4 Proof of Theorem 11.4b = T 3 

In this section, since the approximate solutions in the proof of Theorem 11.31 cannot 
be applied directly to the case of system m (Q in three dimension space, we will 
construct a new approximate system which can be applied to obtain the global weak 
solutions to the three-dimensional system (11.11) (11.31) with h = p and g = 0. 

4.1 A priori estimates 

For constants po and e satisfying 

Po = 50, 0 < e < £\ = min{10 _1 °, 7 / 0 }, 

with 7/0 as in (|1.8I) . we consider the following approximate system 

pt + di v(prt) = evAv + £Ddiv(|Vu| 2 Vu) + ep~ Po , 

< put + pu-Vu- div(pVu) + VP (4.1) 

k = \/idiv(pV7i) + ev\Vv\ 2 Wv ■ Vu — ep~ P0 u — ep\u\ 3 u, 

where v = p 1 / 2 . The initial conditions of the system (HU) are imposed as: 

(p,u)(x,0) = (p 0e ,U0e), (4.2) 

where smooth P-periodic functions po £ > 0 and uq e satisfy 

l|Po£||L 1 ni'v(n) + l|Vpo£ 2 ||L 2 (n) + e ll^ 7 Poe 2 |lL 4 (Q) + e 'll/ J oe ,0 IU 1 (n) — ^ 

and 

/ P 0 e\u 0 e\ 2+VO dx < C, (4.4) 

Jfl 

for some constant C independent of e. 

Let T > 0 be a fixed time and (p, u) be a smooth solution to (14.ip (|4.2I) on P x (0, T], 
Then, we will establish some necessary a priori bounds for ( p,u ). The first one is the 
energy-type inequality. 
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Lemma 4.1 There exists some generic constant C independent of e such that 

r T 


sup [(p\u\ 2 + p + p 7 + sp Po )dx + ( [ p\Vu\ 2 dxdt 

<t<T J Jo J 

f J (|Vu | 4 + |Vu | 2 |«| 2 + |Vy| 4 |?x | 2 + p~ P0 \u\ 2 + p|u| 5 ) dxdt (4.5) 
J J p~ 2po ~ 1 dxdt < C, 


+ £ 


+ £ 


where and throughout this section, for any f, 


fdx = / fdx. 


in 


(4.6) 


Proof. First, integrating (14.11b over yields 

(J pdx)t + e J (| Vu| 2 + |Vu | A )dx = e J p~ Po dx. 

Next, multiplying (I4.1I ) 2 by u, integrating by parts, and using fl4.1D 1 . we have 
i(/ M »*), + J p\V u?i x + VefpWu Pd. + § / p-nu |»<te 

+ £ J p\u\ 5 dx + J u - Vp^dx 

= ^ JvAv\u\ 2 dx + ^ J vdiv(\Vv\ 2 Vv)\u\ 2 dx + e J u|Vu| 2 Vu • Vu • udx (4.7) 

= —^ J\Vv\ 2 \u\ 2 dx — e J vVv -Vu ■ udx — ^ J |Vu| 4 |ri| 2 dx 

< —^ J\Vv\ 2 \u\ 2 dx+ J p\Vu\ 2 dx — ^ J |Vu| 4 |n| 2 dx. 

Then, to estimate the last term on the left hand side of ([4.71) . in a similar way as for 
m , one obtains that for q 7 ^ 1, 


J u ■ Vp q d. 


x = — - 


9-1 

g 

g-1 
1 


9-1 

ge 


g-1 

which, after choosing q = —po, implies 


J p q 4 di v(pu)dx 

J p q ~ l {—pt + evAv + eudiv(|Vu| 2 Vu) + ep~ Po )dx 
( J p q dx) t + q ^~^ £ J p 9 _ 1 |Vu| 2 (l + |Vu| 2 )dx 

p q - l ~ Po dx, 


(4.8) 


6 (po + 1) 


(J p Po dx) t + 


Po(2po + l)e 2 
6 (p 0 + 1) 


J p-^lVul^l + lVul 2 )^ 


+ 

£ 

6 


Po£ 


„-i- 


6 (p 0 + 1) 


2po da 


(4.9) 


< 


Poe 


12(po + 


p Po di vudx 

Yy J p~ l ~ 2po dx + i J p(Vu) 2 dx. 
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Finally, adding (14.61) . (14.71) . and (14.91) together, we obtain (14.51) after using (14.81) . 
Gronwall’s inequality, and the following simple fact that 

^-Po+7-l < p + p-Po. 

Hence, the proof of Lemma 14.11 is finished. 

In the same spirit of the BD entropy estimates due to Bresch-Desjardins 0-13, the 
following estimates also hold. 


Lemma 4.2 There exists some generic constant C independent of e such that 


sup / (|Vu| 2 + e|Vu| 4 ) dx + 
o <t<T J 


(p\Vu\ 2 + p^lVp] 2 ) dxdt 


+ £ 


((Ax) 2 + \S7 v\ 2 \V 2 v\ 2 )dxdt + e 


\S/v\ 4 \V 2 v\ 2 dxdt < C. 


Proof. First, set 

G = evAv + eudiv(|Vu| 2 Vu) + ep ~ P0 , tp e (p) — log p. 
Following the same procedure leading to (12.161) . we can get 

1 + 2 V ^ ( J P~ 1 \Vp\ 2 dx) t + (J u ■ V pdx)t + J p\Vu\ 2 dx 

+ J P\p)p~ l Wp\ 2 dx + {l + y/£) J p _1 G ^Ap-^p _1 |Vp| 2 ^ dx 

= — J Gdivudx + 2 J pVu:Vudx + £ J v\ Vu| 2 Vu • Vu • V log pdx 

— e / p~ Po u ■ V log pdx — £ / p\u\ 3 u ■ V log pdx = ^ I,. 

J J i=i 


Since 

&P- ^P _1 |Vp| 2 = 2vAv, 

the last term on the left-hand side of (14.111) can be rewritten as 

1 


J P G ^A/9- -p |Vp| J dx 
= 2 £ J (Av + div(|Vu| 2 Vu))Audx + 2 £ J p~ Po ~ 1 / 2 Avdx 
= 2 £ J (Av) 2 dx + 2£ J |Vx| 2 |V 2 x| 2 dx + e J |V|Vx| 2 | 2 dx 

+ 2(2p 0 + l)e J p- p °-'\V V \ 2 dx, 

where we have used the following simple fact 

J div(\X7v\ 2 Vv)Avdx = — J |Vu| 2 Vu • VA vdx 

1 


= / |Vu| 2 |V 2 u| 2 dx + 


|V|Vu|Tdx. 


(4.10) 


(4.11) 


(4.12) 
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Then we will estimate each /$(* = !,•■■ , 5) on the righthand side of (14.111) as follows: 


\h\ = 


Gddvudx 


= £ 


< 


v(Av + div(|Vu| 2 Vt’))divztdx + J p po di vudx (4-13) 

f (Av + div(|Vu| 2 Vu)) 2 dx + e 2 J p~ 2po ~ 1 dx + C J p(d\vu) 2 dx , 

|^ 2 1 — /o|Vu| 2 dx + C J p\Vu\ 2 dx , (4-14) 


|/ 3 | = e 

= 2s 

= 2s 
£ 

“ 8 


v|Vv| 2 Vu • V'tt • V log pete 
^ |Vu| 2 Vu • Vit • Vvdx 
dj(\Vv\ 2 djvdiv)uidx 

J \\7v\ 2 \\7 2 v\ 2 dx + Cs J \Vv\ 4 \u\ 2 dx, 


(4.15) 


N =£ 


£ 

PO 


J p Po u ■ V log pdx 
J p~ P0 divudx 


(4.16) 


< Cs 2 / p 2po 2 dx + C p(dxvu) 2 dx, 


\h\=e 


l* 


u\ 3 u ■ V log pdx 


(4.17) 


<CsJ p\u\ 5 dx + Cs J (p + p po )dx + s J \u\ 2 \Vv\ 4 dx. 

Finally, since v satishes 

2vt — sAv — ediv(|Vu| 2 Vu) = —2 u ■ Vu — udivu + sv~ 2p °^ 1 , (4-18) 

multiplying (|4.18l) by s{Av + div(|Vu| 2 Vu) and integrating the resulting equality over 
Q lead to 

s(J |Vu| 2 dx + i J |Vu| 4 dx)t + s 2 J (Av + div(|Vu| 2 Vu)) 2 dx 
(2po + l)e 2 J v~ 2po ~ 2 (\V v\ 2 + |Vu| 4 )dx 

j (Av + div(|Vu| 2 Vt;))udivudx + 2s J (Av + div(| Vv\ 2 Vv))u ■ Vvdx (4.19) 
- j (Av + div(\\7v\ 2 \7v)) 2 dx + C j p(d\vu) 2 dx + ^ j |Au| 2 dx 


+ 


= £ 






















Adding (14.191) to (14.111) . we obtain (14.101) after using Gronwall’s inequality, (I4.12D - 
(|4. 171) . (14.51) . and the following simple fact: 

J (div(|Vu| 2 Vu)) 2 dx 

= f dj(\X7v\ 2 div)di(\X7v\ 2 djv)dx 



The proof of Lemma 14.21 is thus finished. 

With Lemmas 14. II and 14.21 at hand, simliar to Lemma 12.,11 we can prove the following 
Mellet-Vasseur type estimate ( [35]). 

Lemma 4.3 Assume that 7 £ (1,3). Then there exists some generic constant C inde¬ 
pendent of e such that 



(4.20) 


0 <t<T J 

Proof. First, multiplying (14. ID -? by (1 + ln(e + \u\ 2 ))u and integrating lead to 
--y- ( p(e + |u| 2 ) ln(e + \u\ 2 )dx + f pln(e + \u\ 2 )(\Vu\ 2 + ^fe\Vu\ 2 )dx 


2 dt 



(4.21) 


where we have used the following simple facts that 


^ J(e + |u| 2 ) ln(e + |u| 2 )ndiv(| Vn| 2 Vn)dx 



and that 


e 


2 


/ 


p Po (e + |w| 2 ) ln(e + \u\ 2 )dx — e 


1 


p Po (l + ln(e + \u\ 2 ))\u\ 2 dx 



Next, similar to (12.201) and (12.211) . we obtain that 



(4.22) 
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and that 


^ J (e + |u| 2 ) ln(e + \u\ 2 )vAvdx 

< — ^ j |Vu| 2 (e + |u| 2 ) ln(e + \u\ 2 )dx + e J p\Vu\ 2 dx (4.23) 

+ CeJ \Vv\ 2 \u\ 2 dx + ^- J pln(e + \u\ 2 )\'Vu\ 2 dx. 


Finally, it follows from (14.211) (14.231) . (14.51) . (14.101) . and (14.41) that 
/ p(e + \u\ 2 ) ln(e + \u\ 2 )dx < C + [ 

o <t<T. 


sup f p[e + |w| 2 ) ln(e + \u\ 2 )dx < C + f f ln 2 (e + |u| 2 )p 27 1 dxdt. (4.24) 

} <t<TJ Jo J 


Putting (13.51) where a = 1 into (14.241) yields (14.201) . The proof of Lemma B~3l is com¬ 
pleted. 

Next, we will use a De Giorgi-type procedure to obtain the following estimates on 
the lower and upper bounds of the density which are the keys to obtain the global 
existence of strong solutions to the problem (14.11) (14.21) . 

Lemma 4.4 There exists some positive constant C depending on £ such that for all 
(x, t) £ Q x (0, T) 

C- 1 < p(x,t) <C. (4.25) 

Proof. First, it follows from (14.101) . (14.5p . and the Sobolev inequality that 


SUp \\p\\L°° = SUp ||U|| L oo 
0 <t<T 0 <t<T 


< C sup (11^11^2 + || Vu||^4) 2 < C. 

0 <t<T 


(4.26) 


Next, we will use a De Giorgi-type procedure to obtain the lower bound of the density. 
In fact, since w = v~ 4 satisfies 


2wt + 2u ■ Vw — wdivu + ew 2po+3 + 2 ew 1 |V-«;| 2 + 2ew 5 |Vu;| 4 
= eAw + ediv(u; _4 |Vu;| 2 Vu;), 


(4.27) 


_ 1 jc\ 

multiplying Q 4.27 | ) by (w - k)+ with k > ||^(-, 0)|| L oo (f2 ) = || p 0 ' ||L°°(f 2 ) yields that 

sup f (w — k) 2 + dx + e f f(\V(w — k) + \ 2 + w~ 4 \'V(w — k) + \ A )dxdt 
o <t<T J Jo J 

' f ( w\u\\V(w — k) + \dxdt + C ( f (w — k) + \u\\\7w\dxdt (4.28) 

Jo Jn Jo Jn 

C [ [ 1; p~ 4/3 \u\ 4/3 dxdt + ^ [ [ w~ 4 \V(w - k)+\ 4 dxdt, 

Jo Jn k ^ Jo Jn 

where Aj, = {( x,t ) € D x (0, T)\w(x, t) > k}. It follows from Holder’s inequality and 


< C 


< 
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(14.51) that 


J j l Ak p- 4/3 \u\ 4/3 dxdt 

- cr/o 1 ^-) 

c(J J (p + p~ P0 )dxdtj \A k \ 2/3 


4/15 


(4.29) 


< 


< c^\ 


where v k = \A k \. Hence, putting (14.291) into (I4.28[) leads to 

\\(w — fc)+ll L io/3 (nx (o, T )) ^ CvJ , 
where we have used the Sobolev inequality 


(4.30) 


(w — k) 2 + dx + C J J |V(tc — k) + \ 2 dxdt. 

Thus, (I4.30P implies that for h > k, 

i>h < C(h - k)~ w / 3 uf /9 (4.31) 

due to the following simple fact that 

II (w - k)+ |liio/3 (nx( o, T )) >(h- k) 2 \A h \ 3 / b . 

It thus follows from (14.311) and the De Giorgi-type lemma [4(T , Lemma 4.1.1] that there 
exists some positive constant C > C such that 

sup p _1 (x, t) < C, 

(x,t)££lx (0,T) 

which together with (|4.26l) gives (14.251) and finishes the proof of Lemma 14.41 

We still need the following lemma concerning the higher order estimates on (p, u ) 
which are necessary to obtain the global strong solution to the problem (14.11) (|4.2I) . 


(w - A0+lliio/3 {nx( o >T)) < C ™P T 


Lemma 4.5 For any p > 2, there exists some constant C depending on e and p such 
that 

(\\(Pt,Vpt,ut\\ P LP(ti) + ||(p,V^,«)||^ 2iP(n) ) dt < C. (4.32) 

Proof. First, it follows from (14.251) . 04. 5 j) . and (14. 101) that 


sup (|M | L 2 + ||Vu|| L 2 ni 4) + 
0 <t<T 


(|u| 5 + |Vu| 4 |V 2 u| 2 + |Vu| 2 ) dxdt < C. (4.33) 


Next, it follows from (14.181) that v = p 1 / 2 satisfies 

2 vt — ediv((l + |Vu| 2 )Vu) = —div(uu + Vw) — f (u ■ Vv — ev~ 2po ~ l )dx, (4.34) 

|‘ ’| J 
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where for t > 0 , w(-,t ) is the unique solution to the following problem 

J A w = u ■ Vu — ev~ 2po ~ l — p| f(u ■ Vv — ev~ 2po ~ 1 )dx, x € Q, 
W wdx = 0 . 

Since (j4.33[) implies 


u ■ \7vdx 


< ( 'll w llL 2 (a)l|Vu|| Z/ 2( n ) < C, 


we obtain that Vu> satisfies for any p > 2 

l|Vw||z,p(f2) < (7||Au;|| L 3p/(p+3) 

< C(p)\\u\\ LPm \\Vv\\ L 3 m + C{p) 

< C(p)\\u\\ LP (n) +C(p), 

due to (|4.35|) . (I4.33p . and (|4.25[) . 

Setting 

v(x,t) = v(x,t) + 7-^-7 f f (u ■ Vv — ev~ 2po ~ 1 )dxdt, 

^l“l io J 

we get from (14.341) that 

J 2 ht — ediv(|Vh| 2 Vh) = div/, 

= v(x,0), 

with / A e\Jy — uv — Vu>. 

Thus, applying the L p -estimates [H Theorem 1.2] (see also [2][3]) to 
periodic data yields that for any p > 4 


Vv\\% p dt= I ||Vh|| % p dt 


3 P 


JO 

<C{p) ( 1 + f T \\f\\%dt 


'0 


< 


< 


C ^{ 1 + J 0 \\u\\ p LP d?j +C(p) (jf UVull^dt 
C(p) + C(p) 


+ ^J o ll Vt 'll %pdi, 


where we have used (14.361) . (14.371) . (14.251) . and (14.331) . 

Next, note that (ED 2 implies that u satisfies 

u t - + \fe)&u - ^Vdivrt = F. 


where 


F = —u-Vu — p 1 VP + ev 1 // 2 |Vu| 2 Vu • Vrt — ep Po l u — e|ri| 3 u. 


Since 


Ari • Vdi vudx = / |Vdivu| 2 dx, 


(4.35) 

(4.36) 

(4.37) 

(4.38) 

(14.381) with 

(4.39) 

(4.40) 

(4.41) 
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multiplying (14.401) by —2A u and integrating the resulting equality over 9 lead to 


(l|Vii||| 2 )t + /((l + 2 x /e)|Au| 2 + | Vdivn| 2 )dx 
+ 2 ej |u| 3 (|Vu| 2 + 3\X7\u\\ 2 )dx 

= 2 J (it ■ SJu + p~ x VP — eu^ 1 / 2 |Vu| 2 Vu • Vu + £p _P0_1 it) • A udx 

< C\\Au\\ l2 ((||«|| L8 + |||Vu| 3 || i5 )||Vu||^ 5 ||Au|| 3/ 2 5 + ||V«|| La + |M| l2 ) 

< ^l|An|| 2 2 + C(||n||l 5 + |||Vu| 3 ||| 5 )||V«|| 2 2 + C, 

where in the last inequality we have used (|4.33l) . This together with Gronwall’s in¬ 
equality, ()4.33p . and (14.391) gives 

sup ||Vit||| 2 + [ ||V 2 n||| 2 dt < C. (4.42) 

0 <t<T J 

It thus follows from this and the Sobolev inequality that 


IMIl 10 (Ox(o,t)) + ll^ 7 ' u llL 10 /3(ax(o,r)) A C , 

which together with (14.391) (14.42h gives 

\\ u t\\L 2 (Qx(o,T)) + l|V 2 u|| i2 (n x ( 0i r)) + ll-^’[lz,5/2(nx( 0 ,T)) A C. (4.43) 

Using (14.43p and applying the standard L p -estimates to (I4.40p (14.411) MAh with periodic 
data yield that for any p > 2 

ll w t||i,p(ox(o,T)) + l|V 2 -u||z,P(f2x(o,T)) < Cip) + U(p)||P|| iP (nx(o,T))- (4.44) 

In particular, combining M.43I) and (|4.44h shows 

ll^t||_L s / 2 (Qx(0,T)) + H^ 2n llL 5 / 2 (nx(0,T)) - C. 

This combined with M.33p and the Sobolev inequality ( [291 Chapter II (3.15)]) yields 
that for any q > 2 

IMIl8(Ox( 0,T)) + l|Vu|| L 5 (n x ( 0 ,T)) < C(q ), 

which, together with M.39I) and M.41I) . gives 

ll^ ? llL 9 / 2 (ax( 0 ,T)) A C. 

Combining this with (14.441) yields 

ll^tII l 9 / 2 (qx(o,t)) + ll^ 72u lli 9 / 2 (nx(o,r)) - ^ 
which together with the Sobolev inequality ( [291 Chapter II (3.15)]) shows 

IMIz,°°(r 2 x(o,r)) + l|V?x||i45(f2x(o,T)) < C. 

Thus, we get 


H-^llL 4 °(nx(0,T)) < C, 


41 

































which together with (14.441) gives 


IKh||L 40 (nx(o,T)) + I|V 2 «||l40(q x (o i t)) < c. 

The Sobolev inequality ( [29], Chapter II (3.15)]) thus implies 


I|Vu||loo(q x ( 0)T )) < C. 

Then, it holds that for any p > 2, 

\\ u t\\LP(nx(o,T)) + l|V 2 u|| i p(n x ( 0 i T)) < C(p). (4.45) 

With (14.451) at hand, one can deduce easily from (14.181) and (14.21) that for any p > 2, 

ll/ 3 t||-W’(o,T,w 1 >* , (n)) + l|V 2 p||LP(o i 7’ i w 1 .p(n)) < C(p), 

which, together with (14.451) and (14.331) . gives the desired estimate (|4.32l) and finishes 
the proof of Lemma 12.51 


4.2 Compactness results 

We first construct the initial data. Let 


uo = 10 


-10 


(4.46) 


Choose 

o < Poe € C°°(9), \\Vpll 2 \\ 4 L 4<e- 4 "° 

satisfying 

\\P0e - PoWl 1 ^) + \\P0e ~ Po||l7(T2) + II ^ (p)) £ ~ / ? 0^)IIl 2 (O) < e - 
Set 

» = (p§ f+E 2i -) 1/6 . 


It is easy to check that 

Ihn ||poe - polkqa) = 0 (4.47) 

and that there exists some constant C independent of e such that (|4.3D holds. Define 
uo £ as in (I2.58h . It is easy to check that (12.59[) and (14.4[) are still valid. 

Extend then (po £ ,uo £ ) 9-periodically to M 3 . The standard parabolic theory [29], 
together with Lemmas 14.41 and 14.51 thus yields that the problem (14.11) (14.21) . where 
the initial data (po,uo) is replaced by ( po £ ,uo e ), has a unique strong solution ( p £ ,u £ ) 
satisfying 

Pei u ei {Ps)ti V (Pe)t) ( u e)ti V p £ , V p £l V U £ G L P (Q e X (0, T)), 


for any T > 0 and any p > 2. Moreover, all estimates obtained by Lemmas 14.1114.31 still 
hold for ( p £ ,u £ ). 

Letting e —>• 0 + , we will modify the compactness results in Section 12.21 to prove 
that the limit (in some sense) (p, ^fpu ) of (p e , yPp^u £ ) (up to a subsequence) is a weak 
solution to (11.11) (11.31) (11.61) (11.71) with a = 1 and 7 e (1,3). We begin with the following 
strong convergence of p £ . 
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Lemma 4.6 There exists a function p € L°°(0, T; L 1 (0) n L 7 (fl)) such that up to a 
subsequence, 

p e ^ p in L 7 (0 x (0, T)). (4.48) 


Proof. First, for v £ = p^J 2 , it follows from (14.51) and (14.101) that there exists some 
generic positive constant C independent of e such that 


sup f (p £ \u e \ 2 + p £ + pi + ep £ po )dx + [ f p £ \V u £ \ 2 dxdt 
o <t<T J Jo J 


+ £ 


+ £‘ 


(|Vu £ | 2 |« £ | 2 + \Vv £ \ 4 \u £ \ 2 + p e P 0 \u £ \ 2 + Pe\u £ \ 5 ) dxdt 


(4.49) 


J o J pf 2po ~ 1 dxdt < C, 


and that 


sup [(\Vv £ \ 2 + e\Vv £ \ 4 )dx + f f pf 2 \Vp £ \ 2 dxdt 
o <t<T J Jo i 

+ £ J o J ({Av £ ) 2 + \X7v £ \ 2 \V 2 v £ \ 2 + £\X7v E \ 4 \X7 2 v £ \ 2 )dxdt<C. 


(4.50) 


Then, (14.491) and (I4.50jl yield that 


sup ||V / 9 e || L 2 7 /(7+i)(o) < c sup ||p £ /2 || L 2 7(n) sup ||Vp* /2 || L 2 (n) < C, 
0 <t<T 0 <t<T 0 <t<T 


(4.51) 


and that 


e 4/3 [ T \\Vv £ \\%dt < ^ ||Vu £ ||g 3 ( e ||Vu £ || 4 4 ) 1/3 (£||Vu £ ||| 12 )dt 
Jo Jo 

r T 

<C£ \\\\7v £ \\\7 2 v £ \\\ 2 L2 dt<C. 

Jo 

Moreover, note that p £ satisfies 

{p e )t + di v{p £ u £ ) = £v £ Av £ + ediv(u £ |Vu £ | 2 Vu £ ) - e\Wv £ \ 4 + £p~ Po . 
It follows from (14.491) . (14.501) . and (14.521) that 

sup \\p £ u £ \\ L i {n) <C sup \\pl / 2 \\ L i m sup \\p\ / 2 u £ \\ L 2 (n) < C, 

r \S'4-s' r T' 0 C\s'4-s' r T' 


0 <t<T 
rT 


0 <t<T 


' J J {v £ \ Au £ | + v e \ Vu £ | 3 + |Vu £ | 4 )dxdt 

<C£ [ (||« £ || La ||Af, £ || L 2 + K|| L 2 ||Vu £ || 3 6 + ||Vu £ || L 2||Vu £ ||| 6 )dt 
Jo 

< C£ 4 ' 2 + Ce 1 / 3 ^e 4 / 3 J ||Vu £ ||| 6 ( 


< C £ l/3 , 


(4.52) 


(4.53) 


(4.54) 


(4.55) 
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and that 


ef J p £ Po dxdt < e 1 /( 2 P° +1 ) ^g 2 J j p £ 2po 1 dxdt 
< Ce l ^ 2po+1 \ 

The combination of (I4.53[) (|4.56l) implies that 

ll(Pe)t||L 1 (0,T;W / - 1 . 1 (^)) — ^ 

Since 7 < 3, it follows from (14.511) . (14.571) . and the Aubin-Lions lemma that (14.481) holds 
for e —>• 0 + (up to a consequence). The proof of Lemma 14.61 is finished. 

Similar to Lemiria l3.5[ we have the following lemma which deals with the compactness 
of the momentum. 

Lemma 4.7 There exists a function m(x,t) € L 2 (0, T; L 3 / 2 (fi)) such that up to a 
subsequence, 

p £ u £ ^min L 2 (0,T;L p (fl)), (4.58) 

for all p £ [1,3/2). Moreover, 

p £ u £ —>• pu almost everywhere (x,t) € 11 x (0,T), (4.59) 


Po/(2po+l) 

(4.56) 


(4.57) 


where 


u(x,t ) = 


m(x,t)/p(x,t ) 

0 , 


for p(x,t) > 0 , 
for p(x, t ) = 0 . 


Proof. First, it follows from (14.491) . (14.50[) . and the Sobolev inequality that 

j \\^ (P£ U e)\\L 1 (Q,)^' 

<c £ [\\p £ \\ Ll{n) \\pl /2 Vu £ \\ 

L 2 (f2) + \\pl /2 Ue\\ 2 L 2\\Vpl /2 \\ 2 L*(n)) dt 

< C. 


(4.60) 


Next, it holds that 


{p £ u £ ) t + di v(p £ u £ <S> u £ ) - di v(p £ Vu £ ) + VP(p £ ) 

= ev £ A v £ u £ + ediv(u £ |Vu £ | 2 Vu £ <g> u £ ) - e|Vu £ | 4 u £ 
+ v / ediv(p £ Vu £ ) - ep £ \u £ \ 3 u £ . 


For the terms on the left hand side of (14.611) . we have 



PeW 


2 dxdt + 



pjdxdt < C, 



p £ \Vu £ \dxdt < C 



p e \S7 u £ \ 2 dxdt + C 



p £ dxdt < C. 


(4.61) 


(4.62) 

(4.63) 
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Moreover, using (14.491) and (14.501) . we can estimate each term on the righthand side of 
(14.611) as follows: 

£J j (u £ |Au £ ||u £ | + u £ |Vu £ | 3 |'u £ | + |Vu £ | 4 |u £ |)dxcit 
r T 

<Cs \\v £ u £ \\ L 2 (||An £ || L 2 + ||Vu e || 3 6 ) dt 

Jo (4.64) 

+ C (e J J |Vu £ | 4 |u £ | 2 dxdf^ J J |Vu £ | 4 dxdi^ 

< Ce 1 / 6 , 

where in the second inequality we have used (14.521) and (14.551) . 


f/*' 


uA^dxdt 


< Ce 1/5 

< Ce 1 / 5 . 


(17 


1/5 


p £ dxdt J 


■ \ 4/5 

p £ \u £ \ 5 dxdt J 


(4.65) 


Hence, 

||(Pe«e)t||L 1 (0,'7W- 1 ' 1 (Q)) < C. (4.66) 

With (|4.60l) and (I4.66P at hand, one can finish the proof of Lemma 14.71 similarly as 
that of Lemma 12771 

Next, as a consequence of Lemmas 14.61 14.71 and 1 1 .74 similar to Lemma 12.81 one can 
obtain 

Lemma 4.8 Up to a subsequence, 

yfpju £ —> yfpu strongly in L 2 (0, T; L 2 (P)), (4.67) 

with 

yfpu G L°°(0,T;L 2 (P)). (4.68) 

Finally, similar to Lemma f2. 9 1 one can prove the following convergence of the diffusion 
terms. 

Lemma 4.9 Up to a subsequence, 

p £ 'Vu £ -> p a \7u in V , (4.69) 

p“(Vu £ ) tr -> p“(Vu) tr in V, (4.70) 

p“dfvu £ —> p a divu in V'. (4.71) 


4.3 Proof of Theorem 11.41 O = T 3 

On the one hand, for any test function ip, multiplying (I4.53P by ip, integrating the 
resulting equality over P x (0,T), and taking e —>• 0 (up to a subsequence), one can 
verify easily after using (14.481) . (14.671) . (14.471) . (14.551) . and (14.561) that (p,yfpu) satisfies 

unsD. 

On the other hand, let <p be a test function. Multiplying (14.611) by cp, integrating the 
resulting equality over P x (0, T), and taking e -7- 0 (up to a subsequence), by Lemmas 
14.6114.81 and 14.91 we obtain after using (14.63D (I4.65D that (/>, ^pu ) satisfies (11.111) . The 
proof of Theorem 11.41 in the case P = T 3 is completed. 
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5 Proofs of Theorems 11.11 II.3L and 11.4b Cauchy Problem 


Finally, in this section, we indicate how to generalize the approaches in the previous 
two sections to deal with the Cauchy problems in the whole spaces. We start with the 
2-dinrensional case. 

Proof of Theorem li.il - 0 = JR 2 . For do as in (12.5511 and T > 0, we consider the 
system (|1.15l) - (11.17l) in Q e x (0, T) with 

Q £ = (-£- CT V- CT0 ) 2 . 

The initial approximate will be defined as follows. First, choose 
POs e C^(Q £ ) n C' ( f (Q), 0 < po e < e~ 4a ° 


satisfying 


- Po\\lPQ £ ) + \\P0e - Po\\L"t(Q e ) + ||V(po e ^ ~ Po ^llL 2 (Q e ) < £ - 


For v>2 suitably large such that v(a — 1/2) > 5, set 


POe = (p V oi a ~ 1/2) + e 4aou{a - 1/2) ) 


2/(i/(2a-l)) 


(5.1) 


It is easy to check that there exists some positive constant C independent of e such 
that (12.11) with 12 replaced by Q e still holds. Moreover, it holds that 


Jim \\posA ~ PoWmo.) =0, 


(5.2) 


where i)j £ E satisfies 

ipe(x) = 


0, for |x| > 8e CT0 , 
1, for |x| < 4e~ a °. 


Since Pq.viq satisfy (jl.8p . we construct wq £ € Cq°(Q £ ) n Co°(fl) such that 
\\woe - m 0 /p ( 0 1+vo)/(2+r,o) \\ L 2+vo(n) < £ - 


Set 

Then, it holds that 

and that 


1 /( 2 +^ 7 °) 

u 0e — Poe w 0 £■ 


Jim llpoe^Oe - m 0 || L i(n) = 0, 
[ Poe\uoE\ 2+r10 dx < C. 

jQe 


(5.3) 

(5.4) 

(5.5) 


Next, let n = ( 711 , 712 ) denote the unit outward normal to dQ e . We impose the initial 
and boundary conditions on the system (jl.l5j) - (|1.17|) as follows: 

f = 0 on dQe, p{x , 0) = P 0 e{x), X € Q e , 

1 u ■ n = 0, <9iU2 — d' 2 U\ = 0 on dQ £ , u(x , 0) = uq £ (x) x E Q e . 
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It follows from (15.61) 9 that for any smooth function f(x) 


/ fdiUjUiUjdS = / f{diUj — djUi)riiUjdS + / fdjUijiiUjdS 

JdQe JdQe JdQe 


J 

JdQe 


(5.7) 


fuiUjdjUidS = 0 . 


With (|5.7D at hand, for any /, denoting 

[ fdx = f fdx, 

■J J Qe 

one can check step by step that all the apriori estimates in Lemmas I2.1H2.3I and 12.41 
where Q is replaced by Q e , still hold for the smooth solution to the problem (11.151) - 
(11.171) (15.61) . It follows from the L p -theory of parabolic system ( [10]) that Lemma 
12.51 with fl replaced by Q e also holds. Moreover, for T > 0, p € (2,oo), and any 
F € L p (Q e x (0,T)), Theorem 2.1 in [lO] shows that the problem (|2.49l) (15.61) 9 admits 
a unique strong solution u on Q e x [0,T] satisfying 

u,u t y 2 ueL p {Q £ x (0, T)), 

provided p € C(Q e x [0, T]) and p > 0. Therefore, one can use the standard contraction 
mapping principle to prove that for any p > 2 and for suitably small T > 0 the problem 
(|1 . 151) — (11.171) (15.6p has a unique strong solution ( p £ ,u £ ) on Q e x [0,T] satisfying 

Pe, u ei ( Pe)t , V(p £ )t, {u £ )t , V 2 p £ , V 3 p £ , V 2 u £ E L p (Q e x (0, T)). (5.8) 

Then, the apriori estimates obtained by Lemmas 12.41 and 12.51 yield that the local-in¬ 
time strong solution (p £ ,u £ ) in fact is a global one, that is, for any T > 0, the problem 
(11.15j) (11.171) (15.61) admits a unique strong solution ( p £ ,u £ ) on Q £ x [0,T] satisfying 
(15.81) for any p > 2. Moreover, ( p £ ,u £ ) satisfies all the uniform estimates (with respect 
to e) in Lemmas I2.1H2.31 We then extend (p £ ,u £ ) to H x [0, T] by defining 0 outside 
Qe X [0, T], 

Then after routine modifications of the proofs of Lemmas I2.6H2.81 we conclude after 
a standard diagonal procedure that 

Lemma 5.1 There exists a function p € L°°(0, T; L 1 (H) n L 7 (II)) 
subsequence, 

/L -t P in L 7 (0,T;L7 OC (I2)). 

Moreover, there exists some function u(x,t ) such that 

yfpu £ L°°(0, T ; L 2 (II)), (5.10) 

and that up to a subsequence, 

\/~p £ u £ yfpu strongly in L 2 (0, T; Lf oc (0)). (5.11) 


such that up to a 
(5.9) 


Finally, it follows from Lemma 15.II and the proof of Lemma [2.91 that Lemma [2.91 still 
holds for Q = K 2 . With Lemmas 15.II and 12.91 at hand, after modifying slightly the proof 
of the periodic case, one can finish the proof of Theorem 11.11 for the case that D = R 2 . 
The proof of Theorem II.II is therefore completed. 
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We now turn to the 3-dimensional case. 

Proof of Theorem 1 1. /I = M 3 . For cq as i n ([2.551) and T > 0, consider the system 


(11.161) (11.171) in Q £ x (0, T) with 


Qe = (—£ -<T °, £ -<T °) 3 . 


(5.12) 


Define po )£ as in (15.11) . Hence, both (|2.1I) . where Q is replaced by Q e , and (15.21) also 
hold. For a € [3/4,1], set uq £ as in (15.31) . While for a € (1, 2), let 


-1/4 

UQe = P 0e W 0£ , 


(5.13) 


where wq £ € C'^°(H) n Cf^(Q £ ) satisfies 

\\woe-rn 0 /pl /A \\ L 4(Q e ) < £■ 

It is easy to check that (15.41) and (15. 5 p are still valid for uq £ defined in both cases, (15. 3 p 
and (15.131) . Moreover, if a € (1,2), 


' Qe 


POeWoel^dx < C. 


(5.14) 


Next, let n = (ni, 712 , 713 ) denote the unit outward normal to dQ e . We impose the 
initial and boundary conditions on the system (13.21) (11.161) (11.171) as follows: 


/ dn = 0 , u ' n = 0, (V x u) x n = 0 on dQ e , 
\(p,7t)(z,0) = (po £ ,u 0e ), x € Q e . 

Similar to (|5.7D . by (15.15[) i , it holds that 

f djUj'npUjdS = 0 , 


(5.15) 


/ 

JdQe 


(5.16) 


for any smooth function f(x). 

With ([5.16D at hand, denoting that for any /, 


[ fdx= f 
J JQe 


fdx, 


one can check step by step that all the estimates in Lemmas l3.1ll3.4[ where D is replaced 
by Q £ , are still valid for the smooth solution to the problem (13.21) (11.161) (11.171) (j5.15p . 
Therefore, similar as that of two-dimensional case, one can use the standard contraction 
mapping principle and the apriori estimates obtained by Lemmas 13.11 - 13.41 to prove that 
the problem (13.21) (jl.161) (11.171) ([5.151) has a unique solution ( p £ ,u £ ) on Q e x [0, T] for 
any T > 0 satisfying ([5.811 for any p > 2 and all the uniform estimates (with respect 
to e) in Lemmas 13.1113.31 We then extend ( p £ ,u £ ) to SI x [0, T] by defining 0 outside 
Qe x [0,T], 

It then follows from some slight modifications of the proofs of Lemmas 12.61 and 13.51 
13.71 and a standard diagonal procedure that Lemmas 15.11 and 13.81 still hold. With the 
help of these two lemmas, one can adapt the proof of the periodic case to finish the 
proof of Theorem 11.31 for the case that Q = M 3 . The proof of Theorem 11.31 is therefore 
completed. 
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Proof of Theorem HI VL = R 3 . First, we choose the initial approximate as follows. 
For (Jo as in (I4.46H . let 

Qe = (5.17) 

Choose 

0 < Poe £ C^(Q £ ) n c 0 °°(n), ||Vp£ e /2 || 4 4 < e- 4CT0 

satisfying 

IIpoe - PoWl 1 ^) + \\poe ~ Po||z,j(f 2 ) + l|V(po£ 2 _ Pq 2 )\\l 2 {Q) < £ - 
Set 

POe = {Poe + £ 24<j °) 1/6 . (5.18) 

It is easy to check that both (14.31) . where 11 is replaced by Q £1 and (15.21) hold true. We 
then choose uq £ as in (15.31) which satisfies (15.41) and (15.51) . 

Next, let n = ( 711 , 712 , 713 ) denote the unit outward normal to dQ £ . We impose the 
initial and boundary conditions (15.151) on the system (14.11) . Note that (15.161) still holds. 
Moreover, since Vp • n = 0 on dQ e , we have 

g(x)Vp •V(f(x)Vp • n) = 0, on dQ e , (5.19) 

for any smooth functions f{x) and g(x). 

Denoting that for any /, 



one can use (15.161) and (j5.19p to check step by step that all the estimates in Lemmas 
I4.1H4.51 where Q is replaced by Q s , are still valid for the smooth solution to the problem 
(14.11) (15.151) . Then, one can use the standard contraction mapping principle and the 
apriori estimates obtained by Lemmas I4.1H4.5I to prove that the problem (14.1 1) (|5.151) 
has a unique solution (p £ , u £ ) on x [0, L] for any T > 0 satisfying (|5.8I) for any p > 2 
and all the uniform estimates (with respect to e) in Lemmas 14.1114.31 We then extend 
( p £ ,u e ) to D x [0,T] by defining 0 outside Q e x [0,T]. 

It then follows from some slight modifications of the proofs of Lemmas 14.61 and 14.81 
and a standard diagonal procedure that Lemmas 15.11 and 14.91 still hold. With the help 
of these two lemmas, one can adapt the proof of the periodic case to finish the proof of 
Theorem 1 1.41 for the case that D = M 3 . The proof of Theorem ll.4l is therefore completed. 
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